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Abstract 

We consider the transport properties of multiple-particle quantum states in a class of one- 
dimensional systems with a single quantum impurity. In these systems, the local interaction at 
the quantum impurity induces strong and non-trivial correlations between the multi-particles. We 
outline an exact theoretical approach, based upon real-space equations of motion and the Bethe 
ansatz, that allows one to construct the full scattering matrix (S-matrix) for these systems. In 
particular, we emphasize the need for completeness check upon the eigenstates of the S-matrix, 
when these states obtained from Bethe Ansatz are used for describing the scattering properties. 
As a detailed example of our approach, we solve the transport properties of two photons incident 
on a single two-level atom, when the photons are restricted to a one-dimensional system such as a 
photonic crystal waveguide. Our approach predicts a number of novel nonlinear effects involving 
only two photons, including background fluorescence, spatial attraction and repulsion between the 
photons, as well as the emergence of a two-photon bound state. 

PACS numbers: 32.80.-t 03.65.Nk 42.50.-p 72.10.Fk 
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I. INTRODUCTION 



Understanding the properties of a single quantum impurity embedded in a continuum 
of extended states is of central importance in both condensed matter physics and quantum 
optics. In general, a quantum impurity problem is defined by a Hamiltonian of the following 
form: 

H = Hp + Ha + Hmt, (1) 

where Hp describes free propagating quantum particles, Ha describes the internal dynamics 
of the impurity (henceforth we will also use "atom" interchangeably), and i^int describes 
the tunneling processes between the impurity and the free propagating states. In condensed 
matter physics, a notable example of a quantum impurity is described by the Anderson 
Hamiltoniani, where the continuum is the band formed by a free-electron gas, the impurity 
is a local site with a single ci-orbital, and the electrons can tunnel between the impurity and 
the Fermi sea. In quantum optics, the Dicke Hamitonian^, which describes in a full quantized 
fashion the interactions of a two-level atom with photons, also falls into this category. Here 
the extended states are free-propagating photon states, the impurity is the two-level atom, 
and the tunneling term i^int describes the emission and absorption processes. In each case, 
due to the interactions at the localized impurity site, the overall system possesses highly- 
nontrivial and strongly correlated behaviors. 

In this article we focus on the scattering properties of such a quantum impurity, when 
one or more quantum particles are incident upon it. The quantum particles are restricted 
to propagate in a one- dimensional continuum. Such a one- dimensional model is relevant 
to recent experiments on the transport properties of electrons through quantum dots^"^*^, 
and photons through quantum dots^ or trapped atoms^. Moreover, such a one-dimensional 
model can also be relevant for three-dimensional scattering problems. Since the impurity is 
typically far smaller in its spatial extent compared with the wavelengths of incident particles, 
most of the three-dimensional problems involving a single impurity can be mapped into a 
one-dimensional problem, because only S waves need to be taken into account. 

It is known that many quantum impurity problems can not be solved using perturbation 
theory. Instead, since the 1980's, significant efforts have been devoted to non-perturbative 
approaches, such as Bethe ansatz that directly diagonalizes the Hamiltonian^'^'i^iiiii^ii^. 
However, most of these papers assume a periodic boundary condition in order to obtain 
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thermodynamic information of the overall system. Only until very recently were the Bethe 
ansatz approach employed to solve for the scattering properties^i^^. The scattering problems 
involve open boundary conditions and are in fact subtle and require very careful treatments. 

Here we develop a full quantum-mechanically theoretical framework that allows one to ex- 
tract scattering information from the eigenstates of the full interacting Hamiltonian. In par- 
ticular, we emphasize the necessity of the completeness check of the computational scheme, 
in order to obtain the correct description of scattering properties. As an illustration of 
our formalism, the multi-photon problem is completely solved using this approach. The 
formalism, however, is general and can be readily applied to electrons as well. 

The paper is organized as follows: in Sec. [Ill we summarize some of basic results of 
quantum scattering theory. In particular, we discuss the Lippman-Schwinger formalism, 
with emphasis on those aspects that are relevant for our purpose. In Sec. Illll we then discuss 
in details the photon Hamiltonian, and its connections to the Anderson Hamiltonian. Sec. IIVI 
discusses the decomposition of the scattering matrix (S-matrix), which enables and greatly 
simplifies the calculations. Finally, in Sec. |Vl and Sec. IVIIand lVIIl respectively, we present a 
detailed discussion of solving the photon Hamiltonian for its one and two photon transport 
properties. In Sec. IVIIII we briefly discuss the three-particle case. 

II. GENERAL ASPECTS OF THE MANY-BODY QUANTUM IMPURITY SCAT- 
TERING PROBLEM 

Before we begin the mathematical adventures of solving for the quantum impurity scat- 
tering problem, we first provide a brief review of relevant theoretical background. 

A. Revisit of the Lippmann-Schwinger Formalism 

In general, quantum scattering theory deals with a Hamiltonian of the form H = ifo+-f^int, 
where Hq defines the free constituents. Hq, for example, can be Hp + Ha in Sec. HI describing 
particles and the single impurity, ifint defines the interactions between these contituents. 
We will restrict to the case where the interaction range of the quantum impurity is finite in 
space. 

The scattering theory aims to answer the following question: for a given incident multi- 
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particle state, after scattering off tlie quantum impurity, what is the outgoing multi-particle 
state? Long before (t — oo) and after {t +00) the scattering, the incoming and outgoing 
states are away from the quantum impurity and accordingly are outside of the interaction 
range. As a result, both the initial and the final states are free particle states, governed by 
Hq. The quantum impurity therefore can be viewed as an intermediary inducing a mapping 
from one free state to another. The mapping is represented by the S-matrix, which encodes 
all scattering properties. 

In the Lippmann-Schwinger formalism, in order to define the S-matrix, one assumes 
that the interaction Hamiltonian Hi^t was adiabatically "switched on" very slowly from the 
distant past (t — > —00), to its full strength at t = 0, and will be adiabatically "switched off" 
very slowly in the distant future {t —>■ +cxo)^ii^'i^. That is, the interaction i^int is replaced 



The limit e ^ 0+ is to be taken last, after all calculations. The adiabatic switching is 
designed to represent the situation that the incoming particles only interact with the target 
for a short period of time and then fly away^^. 

Let the interacting state of the overall system at t = be The time-evolution of 

the interacting state is described by U(t)|z+), where U(t) is the evolution operator related 
to the time-dependent Hamiltonian Hq + H[^^{t). Following the adiabatic switching on of 
the interaction aforementioned, in the remote past {t —00), U(t)|i^) asymptotically 
approaches e~^^^^^\i) = U°(t)|i), where \i) is a free state. Similarly, with the adiabatic 
switching off of the interaction, in the remote future {t +00), \J{t)\i'^) asymptotically 
approaches U°(t)|/j), where \fi) is a free state. Both \i) and are governed by the free 
Hamiltonian Hq. The subscript i in the state indicates its dependence on \i) and 
directly correspond to the incoming free states prepared, and outgoing free states detected 
in the experiments. Hence they are referred to as "in-" and "out-" state, respectively. 

The three states and as deflned above (Fig. [2]), satisfy the Lippmann- 

Schwinger equations^^ii^: 



by 



HUt)^e-^\'\H, 



e 



(2) 



z+) = + 



1 




(3) 



E- Ho + ie 
1 




(4) 



E- Ho-ie 
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where i , . = is the free retarded Green's function: while ^^-^ — - = is the 
free advanced Green's function. Eqs. ([3]) and (jlj) are apphcable for energy eigenstates, i.e., 
HqIi) = E\i), H\i^) = E\i'^), and Ho\fi) = E\fi). It can be proved that the energies of \i) 
and \fi) are the same as that of 

The S-matrix, in general, is defined as 

S = $^|/.)(^|, (5) 

i 

where the summation is taken over a complete basis of the Hilbert space defined by 
Hq. Once the S-matrix is determined, one can then calculate the scattering properties for 
an arbitrary incident state. For a given in-state |in), the out-state is 

|out) = S|in) = 5^|/,)(z|in), (6) 

i 

and thereby the probability amplitude of finding the out-going particles to be in a state \x) 
is 

(x|out) = (x|S|in) = J](x|/.)(2|in). (7) 

i 

In most practical scattering calculations, one starts with a given \i) and computes 
from Eq. ([3]), and then obtains from Eq. (jlj). By repeating this process for a complete 
set of eigenstates {\i)} of Hq, the S-matrix is constructed. In this route of constructing the 
S-matrix, the unitarity of the S-matrix, i.e., S^S = SS^ = 1, is automatically guaranteed 
by starting with a complete basis set {\i)} for the free Hamiltonian Hq. 

For the impurity scattering problem that we deal with here, however, as it turns out, we 
will be in an unusual situation wherein the interacting state \i~^) is first obtained through 
a non-perturbative Bethe-ansatz technique. In this situation, to construct the S-matrix, 
one has to turn around Eq. (j3j) and (jlj) to compute \i) and \fi) from In doing so, the 
completeness of the set {\i)} thus obtained needs to be explicitly checked, especially since 
the state \i) thus obtained can itself possess rich and entangled structures. (A completeness 
check for typically is far more involved.) This route of constructing the S-matrix, 

and the completeness check, will be explicitly carried out for photon-impurity scattering 
problem in Sec. jVjfor one-photon case, and in Sec. IVII and fVIII for two-photon case. 

As a remark, we note that Eq. (j3j) and (jlj) can also be expressed using the exact Green's 
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funct ionsi^ii^'i^: 

= K) + irV-Hintll) = I'O + G^i^intK), (8) 

H/ — n -\- le 

\^) = m + J, \ . ginti/.) = I/.) + G^i^intl/.). (9) 
£/ — tl — %€ 

Our approach therefore also provides a way to compute the exact Green's functions. 



B. A Simple Example 

As an illustration of the approach of constructing the full S-matrix starting from I?"*"), 
here we give one simple example of a one-dimensional scattering problem wherein a quantum 
particle (or a wave) scatters off a delta potential characterized by Vix) = Vq6{x), as shown 
in Fig. [3l The quantum particle is described by the free Hamiltonian Hq = —(P/dx^. 

An eigenstate of the full Hamiltonian H = Hq + V{x) can be computed straightfor- 
wardly as: 

where = ^'j^^^ , tk = jmy^, and 1 + = 4. To compute \fi), we write Eq. 1^ in the 
real-space representation: 

= + / dx'{x\- ]- -\x')V{x'){x'\i+) 

= {xm + {x\- 1- -|0)\/o(0|z+), (11) 

where (0|i+) = ((0-|i+) + (0+|i+)) /2 = tfc/v^, and Ek = k^. 

Using the fact that the advanced Green's function for Hq = —d'^/dx^ is 

one can easily verify using Eq. flTTl) that the out-state is 

{x\fi) = 'rk%== + tk^^, for all x. (13) 



Similarly, by using the retarded Greens' function for Hq: 



1 la;') = _±_^+ik\x-.'\^ (^4) 



Ek- Ho + ie' ' 2k 
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and Eq. ([3]), the in-state \i) can be obtained as 

{x\i) = , for all x. (15) 

V 27r 

Hence, \i) is a plane wave state \k). Since the set forms a complete basis set of 

eigenstates of Hq: 

Ho\k) = k^\k), 
J2\k){k\ = l, (16) 

k 

the S-matrix is 

s = J2 \f^)('\ = E^'^i - k){k\+h\k){ki (17) 

i k 

with the relations (— A;|S|A;) = r^, and {k\S\k) = tk- 

These relations are consistent with the usual reading of in Eq. ffTOj) . where and 
tk are interpreted as the reflection and transmission amplitude, respectively. The deriva- 
tions here put such an intuitive reading of the interacting eigenstates on a firm theoretical 
foundation. 



C. One- and Two-Particle States Expressed in Second-Quantized Form 

When describing scattering processes involving multiple particles, it is advantageous to 
express the states in second-quantized form. Moreover, similar to the example above, the 
eigenstate in the Bethe ansatz calculations is best expressed in real space representation. 
Here, for convenience, we list the expressions of one- and two-particle states in the second 
quantized form using a real-space representation. 

The basis for real-space representation of one- and two-particle states are 

\x) = c\xm, 

|xi,X2) = -^ct(xi)ct(x2)|0), (18) 
where |0) is vacuum. These states are normalized as 

{x\x') = 6{x — x'), 

{xi, X2\x[, x'2) = l:[S{xi - x[)6{x2 - x'2) ± S{xi - x'2)S{x2 - x[)], (19) 
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with + sign for bosons, and — sign for fermions. ((xi,X2| = 0\-^c{x2)c{xi).) Using this 
basis, any two-particle states |/) is defined as 



with f{xi,X2) = +f{x2,Xi) for bosons, /(xi,X2) = —f{x2,Xi) for fermions. 

f{xi,X2) is in fact the two-particle wavefunction. Let the free Hamiltonian Hq take the 
following form in the second quantization: 



which is relevant to our purpose. One can show that the second-quantized Schrodinger 
equation -ffol/) = E\f) leads to 



the Schrodinger equation in the first quantization {ottof^^. Thus f{xi,X2) has the same 
properties of the two-particle wavefunction when expressed in the first quantization form. 
(Note that similar relations between Eq. fl2T|) and fl22|) hold true, for one-particle wave- 
function, as well as when one-particle external potential, and two-particle interaction are 
includedi^'i^.) 

Below we provide further evidence that f{xi,X2) is indeed a "two-particle" wavefunction. 
For example. 




(20) 




(21) 



Ho{xi) + Ho{x2) f{Xi,X2)=Ef{xi,X2) 



(22) 



{Xi,X2\f) = j dx[dx2-[6{Xi -X[)5{X2 - x'2) ±5{Xi - X2)6{X2 - x[)]f{x[,X2) 
= ^{f{Xl,X2) ± f{x2,Xi)) 
= f{XuX2). 



(23) 



Moreover, for any two-particle states |/) and \g), where 




(24) 



with g{xi,X2) = ±g{x2,Xi) (+ sign for bosons, and — sign for fermions), one has 
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idlf) = I dxidx2dx[dx'2 f* {xi, X2)g{x[, X2)-[S{xi — x[)6{x2 — x'2) ± S{xi — X2)S{X2 — x[)] 
dXidX2^ {f*{xi,X2)gixi,X2) ± f*{xi,X2)gix2,xi)) 
dxidx2f*{xi,X2)g{xi,X2) 

dxidx2{f\xi,X2){xi,X2\g). (25) 

III. THE SYSTEM AND THE HAMILTONIAN 

We now apply the general ideas developed above to the problem of photons scattering off 
a two-level system. Fig. H] shows the schematics of the overall system of interest. The two- 
level system is embedded in a one-dimensional waveguide in which the photons propagate. 
The one-dimensional waveguide can be, for example, a line-defect waveguide in a photonic 
crystal with a complete photonic band gap. A discussion of such a problem for photonic 



crystal experiments is provided in Ref. 15|]. The focus here is on the formalism itself. 



The system is modeled by the Hamiltonian^Si^i: 



d d 
H = I dx{ -iVgc'j^{x)—CR{x) + iVgc\{x)—CL{x) 



+V5{x) ^c|j(a;)(T_ + CR{x)a+ + c\{x)a^ + CL{x)a^ j 

+ E^alae + EgO^^ag (26) 

where Vg is the group velocity of the photons, and c^(x)(c|^(x)) is a bosonic operator creating 
a right-going(left-going) photon at x. V is the coupling constant, a)g{a]^ is the creation 
operator of the ground (excited) state of the atom, 0"+ = alag{a- = ajoe) is the atomic 
raising (lowering) ladder operator satisfying cr+|n, — ) = \n,+) and a^\n,+) = 0, where 
|n, ±) = |n) ® |±) describes the state of the system with n photons and the atom in the 
excited (+) or ground (— ) state. E^. — Eg{= Q) is the transition energy, h is set to 1. 
This Hamiltonian describes the situation where the propagating photons can run in both 
directions, and is referred to as "two-mode" model. 

The aim of this paper is to solve the two-photon transport properties of this Hamiltonian. 
Specifically, we imagine a physical scattering experiment, where two photons incident upon 
a two-level quantum impurity embedded in a one- dimensional waveguide. When the two 



photons arrive at the impurity within a time interval comparable to the spontaneous emission 
lifetime of the impurity, one should expect that the transport properties of the photons are 
strongly correlated, as mediated by the quantum impurity. Here, we develop the theoretical 
formalism to describe such a correlation. 

As a first step, we note that by employing the following transformation 

ct(x)^-i=(4(x)-cl(-x)), (27) 

the original Hamiltonian is transformed into two decoupled "one-mode" Hamiltonians, i.e., 
H = He + Ho, where 



He = J dx{-i)vgcl{x)-^Ce{x) + J dxV6{x) (4(x)cr_ + Ce(x)(j+) + Eea\ae + Egttlag, 

(28a) 

Ho = J dx{-t)vgclix)^Coix), (28b) 

with [He, Ho] = 0. iJo is an interaction-free one-mode Hamiltonian, while He describes a 
non-trivial one-mode interacting model with coupling strength V = \^V. For notational 
simplicity, Vg is set to 1 hereafter. 

As a side note, the Hamiltonian He is closely related to the extensively studied Anderson 
model in condensed matter physics. The Anderson model describes the interaction of the 
conduction electrons with a single quantum impurity- . In real-space, the one-mode Anderson 
Hamiltonian takes the following forn>iii^: 



dx+ / V6{x) \^cl{x)cd,a + cl,^c^{x)j dx + '^ednd^^ + Urid^'ind^i, 

(T (7 

(29) 

where cj.(x) (cct(x)) is the creation (annihilation) operator of the conduction electron with 
spin a, while the operator ^ (Q,cr) creates an electron of spin a on the local impurity at 
a; = 0. rid^a is the number operator of electrons on the impurity. V is the coupling strength. 
ed is the energy of the electron on the impurity, and is degenerate for both spins. 
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The first term in Eq. (129|) describes the kinetic energy of the conduction electrons, and 
the second term describes the interaction ( "hybridization" ) between the conduction electrons 
and the impurity. These two terms closely resemble the first two terms of Hf. in Eq. fl28al) . 
The last term, Una^iria^i, describes the on-site Coulomb interaction between the electrons on 
the impurity. For an isolated impurity, this term has three energy configurations: (i) zero 
occupation with energy Eq = 0; (ii) single occupation by an electron of spin a. The energy 
is = where a =t or I; (iii) double occupation with a spin | and a spin | electron with 
an energy E2 = 2ed + U . When U is large, double occupation is energetically unfavorable. 
In the limit where U +00, double occupation becomes prohibited, and the quantum 
impurity can only accommodate at most one electron, a situation similar to the photon 
Hamiltonian wherein the two-level system can at most absorb one photon at a time. In this 
infinite U limit, since there is no double occupation, the term Und^^Ud^i effectively drops 
out, and the Anderson Hamiltonian Ha is exactly the same as the photon Hamiltonian, H^. 
(except the spin degeneracy). In fact, the general procedures and formalism detailed in this 
article for two photons can be directly applied to the Anderson model for two electrons in 
the spin-singlet state for arbitrary f/^^. Our procedures thus provide a unified computation 
schemes for the transport properties of strongly correlated photons as well as electrons. 

Furthermore, when {Yla^d,a) — 1, ^-c, in the so-called local moment phase, the Ander- 
son model and the Kondo Hamiltonian (s-d model) are equivalent^. Both the Anderson 
model and the Kondo model have recently been applied to nano-structures such as quantum 
dots and single electron transistor-i^>2^>2Ii2^. This connection hints the rich structures of the 
problem of strongly correlated photon transport. On the other hand, unlike the Anderson 
model, where fermionic operators describe electrons, here we have bosonic operators describ- 
ing photons. Consequently, the physics that arises from the existence of a Fermi surface does 
not occur in our system. The transport properties for multi-electrons and for multi-photons 
will be correspondingly different. 

In the following three sections, we will provided a detailed account of our solutions to the 
two-photon transport properties for the Hamiltonian in Eq. (l26l) . 
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IV. RELATIONS BETWEEN THE TWO-MODE AND ONE-MODE S-MATRIX 



The decomposition of the two-mode Hamiltonian H into two decoupled one-mode Hamil- 
tonians [Eq. fl28ap and fl28b|) ] greatly simplifies the calculations. In this section, we present 
a strategy, followed by explicitly detailed calculations, to construct the exact S-matrix of H 
from the scattering properties of He and Ho- 

Since both c^(x) and c|^(a;) can be decomposed into a linear combination of cl{x) and 
cl{x) via Eq. (p7|) . any free one-photon state can be written as 

|^i) = |^)e + |^)o, (30) 

where the subscripts label the subspace spanned by 4(x)|0, — ) or cj,(x)|0, — ), respectively. 
Similarly, since c|j(a;i)c^(a;2), c^(a;i)c^(x2), and c^^{xi)c'^{x2) can all be decomposed into a 
linear combination of cl{xi)cl{x2), cl{xi)cl{x2), and cl{xi)cl{x2), \'^2) can also be written 
as 

1^2) = |^)ee+|^)eo+ 1^)00, (31) 

where the subscripts label the subspace spanned by -^4(a;i)4(x2)|0, — ), cl{xi)cl{x2)\0, —), 
or ^4(xi)4(x2)|0, -), respectively. 

Let the two-mode S-matrix be S, we assert the following decomposition relation: 

S|*l) = Se|^)e + So|^)o, 

S|*2) = See|*)ee + Seo|^)eo + Soo|^)oo, (32) 

where Se is the one-photon S-matrix in the e subspace governed by i^e. So = 1, the identity 
operator, is the one-photon S-matrix in the a subspace governed by Ho- See is the two-photon 
S-matrix in the ee subspace governed by H^, Seo = SeSo, and Soo = 1- Once the terms on 
the right hand side of Eq. ( l32l) are computed, S can be constructed correspondingly. We 
will calculate Se in the next section. Obtaining See|^)ee involves non-trivial calculations, 
and will be done via the Bethe-ansatz approach in Sec. I VII 

To prove the decomposition relation of the two-mode S-matrix [Eq. (132!) ]. we start from 
the asymptotic conditions stated in Sec. Ill Al which relates the in-state \i), out-state 
and the interacting eigenstate \i~^)'- 

U(t)|z+) ) U°(t)|i), (33a) 
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U(t)K+) ) U°(t)|/.), (33b) 

where U''(t) = e"*^"* is the unitary evolution operators for Hq, while U(t) = e~*^* is the 
unitary evolution operator for H. Combining Eqs. (I33al) and (]33bp . one then has 

m = lim ^ lim (v'\tf)V{tf)) (ut(t,)U°(t,)) |z) 

= S|i) (34) 

We note this form of the S-matrix is equivalent to that of Eq. and both have exactly the 
same matrix elements. Also, when \i) is an eigenstate of Hq, Eq. (1341) directly gives rise to 
the Lippmann-Schwinger formalisnt^^. 

To proceed, one recognizes that the photon Hamiltonian H, Eq. fl26|) . can be separated 
a.s H = He + Ho [Eq. (!28l)]. and so is the free Hamiltonian Hq = Hq + Hq, where 

Hq = ldx{-i)Vgc\{x)—Ce{x)+Eealae + Egalag, (35a) 

H° = J dx{-i)vgcl{x)-^Co{x), (35b) 
with [Hq, Hq] = 0. It thus is easily seen that the S-matrix in Eq. ( IMI) can be factored as 

S = ^^hm^^^lim^(uf(t;)Ue(t/)) (ut(t,)U°(t.)) (u°^(t/)U,(t;)) (ut(t,)U°(t.)) 
= SeSo, (36) 

where 

Ue,o(t) = e-'^-°*, 

Koit) ^ e-^^o'°*. (37) 

The factoring of the S-matrix also occurs in situations such as when the Hamiltonian can be 
separated into degrees of freedom of center of mass and relative variables, or the spin and 
spatial coordinates^^. 

The decomposition relation, Eq. fl32|) . follows naturally as a consequence of the factoring 
of the S-matrix. For the one-photon case, we have 

S\i) = E>eSo\i)e + SeSoK)o 

= Se\i)e + So\i)o, (38) 
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where Sg and So take the form as in Eq. ([5]), with \i) restricted to one-photon "e" and "o" 
subspaces. In this derivation, we have used E>e\i)o = \i)o, and E>o\i)e = K)e- Also, for our 
case, E>o = So = 1, as can be seen from Eq. since Ho = Hq. 
For the two-photon case, we have 

S\i) = SeSoK)ee + SeSo|i)eo + §6^0^)00 
= SeK)ee + SeSo|'i)eo + ^(,12)00 

where, again, both See, Seo = SeSo, can be calculated using Eqs. ([3]) - ([5]), and Soo = 1- 

With the decomposition relation established, we now concentrate on constructing Se and 
See for the non-tivial Hg in next two sections. And finally, we will use the solutions of He 
to construct the two-mode two-photon scattering solutions of H = H^, + Ho in Sec. I VIII 



V. ONE PHOTON S-MATRIX: S, 



As a preparation of the two-photon solution, we first briefly summarize the one-photon 
solution for the Hamiltonian H^, which is needed for constructing the two-photon S-matrix 
later. 

One could readily check that the full interacting one-photon eigenstate for H^ takes the 
form22i2i 



dx 



{e{-x) + tke{x))cl{x) 



+ efc(T+ } |0, -) 



= / rfx(/.(a;)4(x)|0,-) +efeCT+|0,-), 



where 



th = 



k-Q-iT/2 



(40) 



(41) 



k-n + iT/2' 

The single photon thus experiences resonance when its energy k is close to the transition 
energy Q of the atom. T = characterizes the width of the resonance and is related 



to the spontaneous emission lifetime of the atom, and 



is the excitation 



^ k~n+ir/2 

amplitude. |0, — ) is the state where there is no photon, and the atom is in the ground state. 
The normalized in-state \k)e and the out-state \ fk)e, constructed from \k~^)e, as shown in 
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Appendix I A 1\ are 

/ rfx05(x)ct(x)|0,-), (42) 



with 

^i{x) = e{x, -\fk)e = -i=e 



ikx 



^%x) = ,{x, -\fk)e = h ( ^e^'" ) = ikH^)^ (43) 



for a// X. The normahzation condition is 

Mk')e = 5{k-k'). (44) 

This demonstrates that one can "read off" the in-state one-photon wavefunction 0i(x), and 
the out-state one-photon wavef unctions, (pfix), respectively, from the "incoming" [x < 0) 
and the "outgoing" (x > 0) part of the photon wavefunction, 0(x). 

Since the set {|A;)e} forms a complete set in the "e" subspace, the one-photon S-matrix 
in the "e" subspace (= Sg) therefore is 

k 

= y^Jk\k)ee{k\ 
k 

= 0'^'^'(^)<^r(^')) [clixM -)(0, -\c,ix')] . (45) 

For any two one-particle states 

= j dxip{x)cl{x)\^,-), 

\x)e = j dxx{x)cl{xM-), (46) 
the transition amplitude is thereby given by 

e(x|Se|V5)e = e{x\k) ee{k\(p) e 

k 

= Y.jj dxdx'hx*{x)<p'^{x)<pf{x'Mx'). (47) 
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We now proceed to solve for the one-photon properties for the two-mode Hamiltonian H. 
Since Ho describes free propagating photons, the one-mode one-photon S-matrix in the "o" 
sector is simply an identity operator in the o subspace, i.e., 

So = l = Y,\^)oo{k\. (48) 

k 

According to the decomposition relation, Eq. (152]) . the two-mode one-photon S-matrix is 

S = Y,tk\k)ee{k\ + Y,\k)oo{k\. (49) 

k k 

For a normalized in-state 

\k')n^ I dx^c^j,{x)\0) 

ik'x 1 r pik'x i 

1 



1 ■k')^ + ^\k')o, (50) 



the out-state is 



S\k')n=^{Se\k\ + So\k')o) 



= ^ \^h\k)ee{k\k% + J2 \k)oo{k\k')o 
= (tk'\k')e + \k')o) 

v2 

= ^(tk> + l)\k')R+^{tk'-l)\-k')L 

= ik'\k')R + fk'\ - k')L, (51) 
where the two-mode transmission amplitude ik' and reflection amplitude ffc/ are 

1 k'-n 

tk' = -ytk' + 1) 



2' ' k'-n + iT/2' 
f^, = -(4, - 1) = (52) 
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respectively, in agreement with previous calculations^Si^i. In the derivations, we have used 

\k')e = J rfx^ct(x)|0,-) 



^ ^ + 10, -) 



\k')o = ^mR-\-k'),), (53) 

and T = V^ = 2V^. 

Fig. [5] plots the transmission and reflection spectrum. On resonance {k = Q), is 0, 
while |ffcp is 1, and the particle is 100% reflected. Note that the effect of the spontaneous 
emission of the two-level system is explicitly included. In the one-dimensional geometry, the 
spontaneous emission therefore does not represent a loss mechanism, and is used here to 
control the coherent transport property of a single photon. 



VI. TWO-PHOTON CASE: CONSTRUCTING See 

Having solved the one-photon case, we now proceed to construct the two-photon S-matrix 
of the two-mode Hamiltonian H. Following the discussions in Eq. ( 132|) in Sec. IIVI the key 
to this is to construct the two-photon S-matrix, See, for H^., which we will undertake in this 
section. 

Before we set out to construct the S-matrix, See, we comment on some of the general 
aspects of two-photon problem that would be useful for this effort. As emphasized before 
in Sec. Ill Al the See matrix is a mapping in the free two-photon Hilbert space. This Hilbert 
space, in its real space representation, consists of all symmetric functions of the coordinates 
of the photons Xi, X2, and is spanned by a complete basis {\Sk^p)ee : k < p} deflned as 

\Sk,p)ee = jj dXidX2Sk,piXi,X2)-^cl{Xi)cl{x2)\^,-), (54) 



with 



{Xl,X2\ X2 



1 1 



ikxi ipx2 I ikx2 ipx 



e 



_|_ Qikx2^ipxi'^ 



27rv^ 

e'^"^ cos (Ax) , (55) 



27r 
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where 

E^k + p (56) 

is the total energy of the photon pair, 

Xc = {Xi +X2)/2, 

X = xi — X2, (57) 

are the center of mass coordinate and the relative coordinate, respectively. 

A = {k-E/2) = {k-p) /2, (58) 

measures the energy difference between two photons. The completeness of {\Sk,p)ee '■ k < p} 
is expressed by 

/OO fp 
dp I dk\Si;^p) ee ee{Sk,p\ 
-OO J —OO 

2 poo r-\-oo 

= -/ dp dk\Sk,p)eeee {Sk,p\, (59) 

^ J -OO J-oo 

using the fact that \Sp^k) — \Sk,p)- As a side note, in computations related to the two-photon 
Hilbert space below, we will adopt two equivalent set of variables: {k,p), with k < p, and 
(A, E), with A = {k — p)/2 and E = k+p, and use the two sets of variables interchangeably. 
The Jacobian between the two sets is 1. Thus, 

p poo poo pO 



dp dk= dE dk= dE dA. (60) 

— OO J —OO J —OO J —OO J —OO J —OO 

Alternatively, the same Hilbert space can instead be spanned by another basis {|y4fc,p)ee '■ 
k < p} defined as 

\Ak,p)ee = JJ dxidx2Ak,p{xi,X2)^cl{xi)cl{x2)\0,-), (61) 



with 



,{X1, X2\Ak,p)ee = Ak,p{xi, X2) 



^-^sgn(a;) (e^^^^^e^^'^^ - e^'^^^e^f^i) 

277 



sgn(x) e'^""' sin (Ax) , (62) 
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where sgn(x) = 9{x) — 9{—x) is the sign function. By definition, one has \Ap^k)ee = — |^fc,p)ee- 
We emphasize that, while both {\Sk,p)ee : k < p} and {\Ak^p)ee : k < p} are complete^, 
arbitrary hnear combination {ak^p\Sk,p)ee + ^A:,p|^/c,p)ee '■ k < p} may not be. The properties 
of the two complete sets {\Sk,p)ee '■ k < p} and {|v4jtp)ee : k < p} are summarized in 
Appendix [Bl Here we only emphasize that the two states |S'fc^^p^)ee and \Ai:2^p.^)ee are not 
orthogonal to each other. 

We now proceed to construct See as follows: we start in Sec. IVI Al by deriving the real- 
space equations of motion from the Schrodinger equation H^l^) = E\^). In Sec. IVI Bt we 
then solve the real-space equations of motion using the standard Bethe-ansatz approach to 
obtain a class of eigenstates of the interacting Hamiltonian Hg. In Sec. IVI CI we obtain the 
corresponding "in-" and "out-" states from the interacting eigenstates using the Lippmann- 
Schwinger formalism discussed in Sec. Ill Al Through a completeness check, we show that the 
in-states thus obtained are in fact not complete. Instead, in order to span the two-photon 
free Hilbert space, one must supplement it with another class of states: a two-photon bound 
state. In Sec. IVI D I we show that the two-photon bound state is also an eigenstate of See- 
This, together with the completeness check upon the two classes of solutions, allow us to 
determine the exact form of See- Finally, in Sec. IVI Et we discuss some of the properties of 

See- 



A. Equations of motion and boundary conditions in real space 

An eigenstate for H^, has the general form: 

= 1$) = dxidx2g{xi,X2)-^cl{xi)cl{x2) + j dxe{x)c\{x)a^ |0, (63) 

where e{x) is the probability amplitude distribution of one-photon while the atom in the 

excited state. Due to the boson statistics, the wavefunction satisfies g{xi^X2) = +g{x2,xi), 

and is continuous on the line Xi = X2- 

From He\^) = E\^), by equating the coefficients of 44|0, — ) and 4cr_|_|0, — ), respectively, 

we obtain the equations of motion: 

d d \ V 
'^dx ^'dx ^ ) "71 (e(a;i)(5(x2) + e{x2)S{xi)) = 0, (64a) 

-t^ -{E- n)^ e{x) + ^ (^?(0, x) + g{x, 0)) = 0, (64b) 
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where g{Q,x) = 1/2 x igiO-,x) + giO+,x)) = gix,0) = 1/2 x {g{x,0-) + gix,0+)). 
The functions g{xi,X2) and e(x) are piecewise continuous. For any such piecewise con- 
tinuous function f{x), the derivative of f{x) is the ordinary derivative plus contribu- 
tions from the jump discontinuities. If Xq is such a discontinuity, we add a term^- 
hm^_^+ f{x) - hm^_^- /(x)J 5{x - Xq). 
The interactions occur on the coordinate axes: xi = 0, and X2 = 0. Applying the 
equations of motions on the boundaries between adjacent quadrants gives the following 
boundary conditions on the boundary of quadrants II and III (xi < 0): 

- i {g{xi, 0+) - g{x,, 0-)) + ^e(xi) = 0, (65a) 

(^-^^ -{E- n)^ e(xi) + ^{gi^i, 0+) + g{xi, o-)) = O, (65b) 
and on the boundary of quadrants II and I (x2 > 0): 

- t (^?(0+, X2) - g{0-, X2)) + ^<^2) = 0, (66a) 

-z— -{E- n)\ e{x2) + ^(^(0+, X2) + ^7(0", X2)) = 0. (66b) 
These boundary conditions must be supplemented by a further condition 

e(O-) = e(0+), (67) 

which arises directly from Eq. fl64bp and ensures the self-consistency. When there are more 
than two photons, this condition gives rise to the Yang-Baxter relationiii^. 

The xi-axis, X2-axis and the line xi = X2 dissect the X1-X2 plane into six regions (Fig. [6l). 
When g{xi, X2) is given in either one of the six regions, one could use the boundary conditions 
to obtain g{xi,X2) in all other regions. 

For example, at the boundary between quadrant II and III, using Eq. fl65ap . we have 

g{xi, 0+) = g{x^, 0-) - ^^e{x^). (68) 



Substitute this into Eq. ( I65b[) . one obtains 



- ^^e(xi) = (^{E -n)+ 20 e(xi) - V2Vgix,, 0"), (69) 



which has the solution 



e(xi) = ce+<(^-^)+*^)^^ +2e+<(^-^)+'^)"i T e-<(^-^)+^?)-' f-V2Vg{x,0-)) dx, 

J —00 



(70) 
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where c is an integration constant. By requiring e{x) to be zero when the couphng strength 
V is zero, we have c = 0. Hence 

e{x,) = ^e+<(^-^)+*^)-^ r e-<(^-f^)+'5)- (-V2Vgix, Q-)) dx. (71) 

J — oo 

B. Constructing eigenstates of He using Bethe ansatz 

We now solve the two-photon equations of motions of in Sec. IVI Al using the Bethe 
ansatz. The Bethe ansatz usually postulates that the eigenstates are superpositions of a few 
extended plane waves when all particles are away from the impurity^^i^. We shall also call 
these solutions of the Wiegmann- Andrei states, after the two authors who worked out 
similar solutions for the Kondo model^*^ and for the Anderson modeU^. 

For the two-particle case, the Bethe ansatz postulates that in regions 1, 2, and 3 (Fig. Ej), 
the two-photon wavefunction has the form 

Q^^ikx^+ipx2 ^ j^^^ipx,+ikx2^ in region 3 (xi < < 0); 

g{xi, X2)= < 52e*^^'i+*P^2 + j^^^ipx^+ikx2^ in ^ggjQ^ 2 (xi < 0, Xa > 0); (72) 

Q^^ikx^+ipx2 ^ j^^^ipx,+ikx2^ in region 1 (xs > Xi > 0). 

The wavefunction in other regions is defined by boson symmetry. The goal of the compu- 
tations, based upon the Bethe ansatz, is then to check that such a form indeed satisfies the 
appropriate equations of motion, and in the process of checking, to determine all constraints 
relating the A's and -B's coefficients. Since by construction, g{xi,X2) already satisfies the 
equations of motion in regions 1,2, and 3: 

for xi 7^ and X2 ^ 0, all we need is to use the boundary conditions [Eq. (l65l) and fl66l) ] and 

the self-consistency condition [Eq. fl67|) ] to determine the constraints on A's and -B's. 

At the boundary between quadrant II and III, since in region 3 (xi < X2 < 0), g{xi, X2) = 
^^^ikxi+ipx2 _^ ^^^ipxi+ikx2 ^ ^l^gj^ l^as 

^(xi, 0-) = Ege*'^^! + Ase'P''^ (74) 
Therefore, using Eq. (!7T|) . we have, for x < 0, 

21 



Plugging e(x) to Eq. (l68i) . we obtain 

_ V 
5,(xi,0+) = 5((xi,0 )-i-j^e{xi) 

- + V2 ^ \v-^ + ^ k-n + zr/2 J 

= tpS3e^'^^^+4A3e^P"\ (76) 

and therefore, in the whole quadrant II [xi < 0, 0:2 > 0), using the Bethe ansatz form of 
g{xi,X2) [Eq. (!72l) ]. we have 

B2 = tpB^] A2 = tkA^. (77) 

One can understand this expression by realizing that when going from quadrant III to 
quadrant II, Xi is unchanged, while X2 : 0^ ^ 0+. Consequently the part of the wave 
function B^e^^^^^"^^^"^ acquires a transmission coefficient tp, and the part of the wave function 
j^^Qipxi+ikx2 acquires a transmission coefficient tk- 

In addition, from the expression of e(x < 0) [Eq. (175|) ]. we have 

We apply the same procedures to the next boundary. The boundary conditions on the 
boundary of quadrants II and I (x2 > 0) are (reproduced here from Eq. (!66!) ): 

-I (^(0+, X2) - (7(0-, 0:2)) + ^e(a;2) = 0, 

-I— -{E- ^]) j e(x2) + ^(^?(0+, 2:2) + ^?(0-, X2)) = 0. 

As previously, from the first equation, we have 

^(0+, X2) = ^7(0", X2) - ^^e(x2). (79) 

Substitute into the second equation, we obtain 

- ^^e(x2) = (^E «0 e(x2) - 72^^(0", X2). (80) 

Since 

^(0-,X2) = tpEse*^^^ +4A3e^'="S (81) 
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we have, for x > 



e[x) -y^V + p - n + ^^/2 J ' ^^^^ 



and 



_ V 
^(0+,X2)=^(0 ,X2)-i-j=e{x2) 



(,,B3e- . «.4,e-.) + (-,r) (^^^ + ^ _ ^ ^ 



1 - 7^ r- 1 t^B^e'P''^ + ( 1 tkB-^e'^''^ 

k-n + iT/2) ^ ^ \ p-n + iT/2) ^ 

= t,h[B,e'^^^ + A,e'^^^). (83) 

Therefore, in region I (x2 > xi region of quadrant I), using the Bethe ansatz again, 

Bi = tptkBs; Ai = tptkAs. (84) 

One can understand this expression by reahzing that when going from quadrant II to quad- 
rant I, X2 is unchanged, while Xi : 0^ — 0+. Consequently the part of the wave func- 
tion tpi?3e*^^^"*"*^^'2 acquires a transmission coefficient tk, and the part of the wave function 
tfcAse*^^^^*'^^^ acquires a transmission coefficient tp. 

Also, from the expression of e(x > 0) [Eq. flS^ ]. one has 

Combining Eqs. fl75]) and flS^ . together with the self-consistency condition e(0~) = 6(0"*") 
[Eq. (l67j) ]. we can determine the ratio of B3/A3 from 



^p-Q + iT/2 k-VL + iV/2J \k-VL + iVl2 p-Q + iT/2j' 
which simplifies to 



(86) 



B^ k — p ~ iT 

A's ^ k-p + iV ^ 
As can be seen from Eq. ( 163!) . the two-photon wavefunction g{xi, X2) and the amplitude e{x) 
completely determine the interacting eigenstate of H^. Fig. [7] summarizes the two-photon 
wavefunction g{xi,X2) in the entire X1-X2 plane, as well as e{x) for all x. 
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We now need to extract the information of the in- and out-states from the interacting 
eigenstate As shown in Appendix [K\ the in-state \i) and the out-state are 

\i) = j dxidx2gi{xi,X2)^cl{xi)c\{x2)\^,-), 

\fi) = j dxidx2gfixi,X2)-^cl{xi)cl{x2)\^,-), (88) 

where, for all Xi and X2, 

gi{xi,X2) = g{xi <0,X2 < 0) [g{xi,X2) in quadrant III], 

gf{xi,X2) = g{xi > 0,X2> 0) [g{xi,X2) in quadrant I]. (89) 

Note this result is consistent with the intuitive notion that the in-state is the "incoming" 
part, i.e., the Xi,X2 < region, of the full interacting state; while the out-state is the 
"outgoing" part, i.e., the Xi,X2 > region, of the full interacting state. 

The in- and out-states, when explicitly spelled out in real-space, have non-trivial struc- 
tures. g{xi,X2) in the full quadrant III {xi,X2 < 0) is, 

g{xi,X2) = (Age^'^^i+'P^^ + B3e'P^'+'''^^) 9{xi - X2) + (^ge^'^^^+^P^^ + A3e'P'''+'''''^) 9{x2 - xi) 
(x[{k-p + zr)e^^^i+*P^2 + {k~p- zr)e**'^^+^'=^2] 9{xi - X2) 
+ [{k-p- ir)e^'=^^+^P^^ + {k~p + ir)e^P^i+*'=^2] e{x2 - xi) 

= {k-p) (e^^^'i+^P^2 ^ ^^kx2+^px,^^ ^ (^^ikx,+ipx2 _ ^ikx2+ipx^^ ^^^(^^^ _ 

(x{k~p) Sk,p{xi, X2) + iTAk,p{xi,X2), (90) 

where Sk,p{xi,X2) and Afcp(xi,X2) are defined in Eq. (155|) and (162!) . respectively. Therefore, 
in the entire X1-X2 plane, the in-state photon wavefunction gi{xi,X2) is 

gi{xi,X2) (x{k-p) Sk,p{xi,X2) + iV Afc,p(xi, Xa). (91) 

Similarly, in the entire X1-X2 plane, the out-state photon wavefunction gf{xi,X2) is 

gf{xi,X2) oc tktp [{k - p) Sk,pixi,X2) + iT Ak.p{xi,X2)] . (92) 

Note that g{xi,X2) is equal to zero when k = p in the entire X1-X2 plane. 

C. In- and out-states from the Wiegmann- Andrei state 

Following the discussions of the previous section, we therefore define 

\Wk,p)ee = {k - p)\Sk,p)ee + i^\Ak,p)ee, (93) 
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and discuss some of the general properties of \Wk,p)ee- These states are obviously impor- 
tant for the scattering problems, since they are the eigenstates of the S-matrix, See, with 
eigenvalues tktp, as can be seen from Eqs. (I9T1) and ( l92l) . Each state therefore is directly 
analogous to a so-called "scattering channel" in the partial wave expansion^^'^^. Below we 
will normalize these states and show that they are orthogonal to each other (as expected, 
since they are, after all, eigenstates of the S-matrix with different eigenvalues). Most impor- 
tantly, and perhaps surprisingly, even though they directly arise from the standard Bethe 
ansatz approach, they are in fact incomplete and thereby can not span the free two-photon 
Hilbert space. 

From the definition, Eq. ( l93l) . it is clear that 



\Wk,p)ee = 0, when k = p, (94a) 

\Wp^k)ee = -\Wk,p)ee, foT anj k, p. (94b) 
The normalization and the check for orthogonality is straightforward: 

ee{^^ki ,pi I ^^k2,P2 ) ee 
= (fci — pi) {k2 ~ P2) ee{Ski,pi \Sk2,P2)ee ~^ ^ ■ ee(^fci,pi \^k2,P2)ee 

+ (^1 ~ Pl)(^r) ■ ee{Ski,pi\Ak2,P2)ee + (^2 ~ P2){^i^) ' ee{^ki,pi\Sk2,p2) ee 

= {ki -pi){k2 -P2) [5{ki - k2)6{pi -P2) + 5{ki - p2)S{pi - ^2)] 

+ [6{ki - k2)6{pi - P2) - 6{ki - p2)6{p, - k2)] 
= [{k, - pi){k2 - P2) + 5{k, - k2)5{p, - P2) 

+ [{ki -p,){k2-p2) - r^] 6{ki-p2)6{p, - k2), (95) 

where we have used the overlap between various |S')ee and |v4)ee states, as provided in 
Appendix [Bl One thus is led to the definition of \Wk,p)ee'- 

l^'^.^)- - -;y^^=r^\Wk,p)ee = 7|^i^l^'=>.)ee 

^ i2A\Sk,p)ee + lT\Ak,p)ee), (96) 



V4A2 + r2 

with \Wk,p)ee being normalized to 



'.{Wk^,pi\Wk2,p2)ee = 5{ki - k2)5{pi-p2), (97) 
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when ki < pi and ^2 < P2, or ki > pi and ^2 > P2- For other cases, one has 

ee{Wk,,p^\Wk^^p^)ee = -S{ki - k2)5{pi - P2) , (98) 



which arises from Eq. fl94bp 



Using the normahzed \Wk,p), the in-state and out-state thus are 

\i) = \Wk,p)ee = J dXidX2Wk,p{Xi, X2)^cl{Xi)cl{x2)\(/i, -) , 
\fi) = tktp\Wk,p)ee = j dXidX2tktpWk^p{Xi,X2)^cl{Xi)cl{x2)\^,-), (99) 



with 



Wk,p{Xi,X2) = ee{Xl,X2\Wk,p)ee 



[{k - p) Sk,p{Xl,X2) + iV Ak,p{Xi,X2)\ 



^{k-pf + T^ 

e'^""^ [2A cos(Ax) - rsgn(x) sin(Ax)] . (100) 



2n 

In describing the scattering process, we will need to find all the eigenvalues of the S- 
matrix. The set of these eigenstates then span the free two-photon Hilbert space. To 
check whether {\Wk',p')ee '■ k' < p'} is complete, one could start with an arbitrary state, for 
example, \Sk,p)ee, project out all \Wk'y) components and calculate 

\^k,p) = \Sk,p)ee — ee{Wk',p'\Sk,p)ee \Wk',p')ee- (lOl) 

k'<p' 

If the set {\Wk',p')ee '■ k' < p'} were complete, such a computation should yield \Sk,p) = for 
arbitrary k and p. This computation is performed in Appendix O Surprisingly, independent 
of the choice of \Sk^p)ee, the computation results in \Sk,p) oc \BE)ee, where 

\BE)ee = j dx,dX2BEix,,X2)^cl{x,)clix2)\fD,-), (102) 

with 

/r 

ee{xuX2\BE)ee = ^^(Xl, X2) = ^e^^^^^l^l, (103) 

v47r 



and normalized as 

,,{BE'\BE)ee = 6iE-E'), (104) 

where Xc = {xi + X2)/2, and x = xi — X2- The defining feature of Be{xc,x) is that, when 
X ±00, \Be{xc,x) \ 0, and therefore \Be) is a two-photon bound state. 
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The set {|W/c,p) : VA; < p}U{|-^-b)} together forms a complete basis of states, and any 
symmetric functions of xi and X2 can be expanded using {Wk,p{xi, X2), Be{xi, X2)}- The 
completeness of this basis is crucial for discussing the transport properties of scattering 
problems. 

D. Two-Photon bound state is an eigenstate of the S-matrix 

We now show that the two-photon bound state is an eigenstate of the S-matrix, with 
eigenvalue 

E-2Q-2ir , , 

= g-2n + 2.T - 

This therefore concludes the calculations of the S-matrix. 

Suppose that in region 3 {xi < X2 <0), g{xi,X2) takes the following form 

g{xi, X2) = = gi(fi-ir)a.i/2gi(E+ir)x2/2_ ^^Qg^ 

We then apply the same procedures as previously to obtain g{xi,X2) in any other regions. 
One first has 

g{xi < 0, 0") = e'^^i/'+^^i/^ (107) 
With the same boundary conditions between quadrant III and quadrant II, we have 

e{x, < 0) = ,e+*(^-"+'^/')^i r e~'^''-''+'''/''^^'{-V2V)e'''^'/^+^^'/^dx' 

J —00 

_ 2V2V iExi/2+Txi/2 (.r.n^ 

~E-2fl + 2ir ' ^ ^ 

where we have used V"^ — F. Note that the resonance occurs at £■ = 2Q. This is in contrast 
to the single particle excitation in {|Wjfc,p)ee} where the resonances occur at k — or p — D,. 
Proceed as before, 

g{xi, 0+) = -i^e{xi) + g{xi, 0~) 

^ ■ V 2\/2V tExi/2+ra:i/2 , iExi/2+ra;i/2 

^E-2Q.^2iT 

E-2n ^iExi/2+rxi/2_ ^-^Qg-) 



E-2n + 2ir 

Therefore, in quadrant II {xi < < X2), in accord with the Bethe ansatz, we postulate 
g{x, < 0, X2 > 0) - ^-2^ gigxe+r./2_ (^^q) 
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To extend to quadrant I, we first obtain g{0 ,X2>0): 



thus 



^(0-,x2 > 0) = _^^^e'E.,/2-r.,/2 

' ' E-2Q + 2iT ' ^ ' 



e{x2 > 0) 

+i{E-n+ir/2)x2 ^-~i{E^Q.+iV I2)x' ^ \^V) ~ ^lEx' /2-Vx' /2^^l 

7_oo E — 2VL + 2iV 

^iEx2/2-T/2x2_ (^112) 



le 



2^/2V ^iEx2/2-T/2x2 



E~2n + 2iT 
From tliis, we obtain 



5f(0+,X2 > 0) = -2-^e(x2) +5'(0 ,X2) 



.^^I ^ E -2^1 \ ^iEx2/2-rx2/2 



E-2n + 2iT E-2n + 2iT 
E -2^l-2iT ^■^^^/2-rx2/2^ ^^^^^ 



E-2Q + 2iT 

and thus, in region I (x2 > Xi > 0), applying the Bethe ansatz again, 

Therefore, in the full quadrant I, 

gix^ > 0,X2 > 0) = ^ - 2^ - ,^,^^r./2g(_^^ E -2n-2zT ^^^^_^^/ 

' ^ E-2n + 2iT ^ ^ E-2n + 2iT ^ ^ 

^ E — 20, — 2ir ^-^^^_p|^|/2 

E - 2fi + 2ir 

= t£;e'^"^-^l"l/l (115) 

While in the full quadrant III, 

g{xi > 0,X2 > 0) = e^^^^+^^/2^(-x) + e*^^=-^^/2^(x) 

= (116) 

Finally, note that for the two-photon bound state, the self-consistency condition e(0~) = 
6(0"*") is automatically satisfied. This proves that See|-B_E)ee = tE\BE)ee, and therefore \BE)ee 
is an eigenstate of See- Fig- [8] summarizes g{xi, X2) in the entire X1-X2 plane, as well as e(x) 
for all X, for the two-photon bound state \Be)- 
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E. The S-Matrix for S 



From the definition of tlie S-matrix, See = ^|in) |out)(in|, tlie two-photon one-mode 
S-matrix therefore is 

ee ee {Wk,p\+Y.^E\BE)eeee{BE\. (117) 

fc^p E 

For |in) = \Wk,p)ee, or |Ss)ee, the out-state |out) = tktp\Wk,p)ee, and tE\BE)ee, respectively. 

It should be explicitly pointed out that the S-matrix defined above, Eq. flll7p . describes 
the physical scattering process that the photon in-state is mapped to the out-state via 
|out) = See|in). This definition of the S-matrix is exactly the same as that in the usual 
scattering theory. In the literatures on Bethe ansatz, unfortunately, sometimes a different 
definition is adopted^'^. There, the S-matrix is defined to be Eq. (187|) . the relative phase 
of the two plane waves of the wavefunction in region 3. 

Below we summarize several computations that are needed for two-mode calculations 
later. The details for these computations are provided in Appendix [Dl We first mention the 
results of ee{Sk2,p2\^ee\Ski,pi)ee, the momcutum distribution of of the out-state ee(5'fc2,p2 l^^t) 
for in-state \Sk-^^p-^)ee in ee subspace: 

ee(5'fc2,p2|See|5'fc,,p,)ee = tk^tpJih - k2)6{pi - P2) + tkj.p^5{ki - P2)S{k2 - Pi) + B6{Ei - E2) , 

(118) 

where the first two terms of product of delta functions indicate the uncorrelated part of 
the S-matrix, which are simply the direct and exchange terms of each individual incident 
momentum, and can also be written as ^^^^^^^(Ai — A2)5(-E'i — -E2) + tkitpiS{Ai + A2)S{Ei — 
E2). The third term 

16ir2 Ei-2Q + ir , , 

B = fll9) 

71 [4Al-{Ei-2Q + try][4Al-{Ei-2Q + ir)^y ^ ' 

in contrast, indicates the strong correlations between the two photons, and manifests as 
the background fluorescence due to the scattering. Note that this term does not conserve 
individual energy of each photon, but only the total energy. When Ai 7^ A2, \B{Ei, Ai, A2) P 
is the probability density for the outgoing photon pair in (E'l, A2) state, when the incoming 
photon pair is in (ii^i, Ai) state. 

The uncorrelated part in Eq. (11181) comes entirely from the first term in Eq. (11171) . 
Ylik<p^k'tp\Wk^p) ee ee{Wk,p\] whilc the Correlated part in Eq. (I118p . B6{Ei — E2), has contri- 
butions from both |Vrfe,p) and \Be) in Eq. (I117p . 
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For the same in-state |m) = \Skx^p^)ee = \SEi,Ai)eey could also write down the real- 
space representation of the out-state: 

ee('2^C7 *^|out)ee ee (*^C7 -^l See | i-^) 

= ^ ] SE2,A2i-'^c, x)ee{SE2A2\^ee\^Ei,Ai)ee 
E2,A2<0 

= Se2,A2{^c, x) {tk^tpx ee{S E2,A2\^ E-i_,Ai) ee + B6{E2 — Ei)) 

E2A2<0 

= tkj^tpj^SEi,Ali-'^cy ^) + BSEi,A2i.^cy ^) 

A2<0 

_ iEiXc^ ( . . cn^(\,r\ — i{Ex-2Q)\x\/2~T\x\/2\ 

- e V ^^^^^ 4Af - [E, -2n + iTY^ ) 

(120) 

which takes the form e*'^^^'=(x|0), where is the wavefunction in the relative coordinate 

X. The deviation of the out-state wavefunctions from that of interaction-free case is large 
when Ai ~ ±(i?i/2 — fi), i.e., when at least one of the incident photons is close to resonance. 

VII. TWO-PHOTON CASE II : TWO-MODE MODEL 

We now compute the two-mode two-photon scattering properties. To analyze a two- 
photon scattering experiment, one first projects the wave packets describing the two photons 
to each IS'^p), and applies the previous discussions to each component. Specifically, consider 
an in-state 

|in) ^ \Su,,p:)rr = J dx.dx^:^ (e^fci-i+^PiX2 + g.fc.x.+.p...) _Lct^(a;04(x2) |0, -), 

(121) 

which describes two incident photons of plane waves from the left with momenta k and p 
respectively. To apply the decomposition relation, Eq. (132!) . we first decompose the in-state 
|'S'fci,pi) to the components in ee, oo, and eo subspaces, followed by computing the scattering 
states in each subspace, and finally transform the results back to RR, LL, and RL spaces. 
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The two-mode out-state thus obtained is (please refer to Appendix [E] for details) 



|out) = S|in) 



dxidx2 t2(a;i, 0:2) ^c|j(a;i)cy 0:2)10, -) 



+ j dxidx2 r2{xi,X2)-^c{{xi)c[{x2)\^,-) 

+ / dxidx2rt{xi,X2)c^n{xi)c[{x2)\(l>,-) (122) 



where 



t2[XuX2)-e yt,,tp, cos [IX,X) 4^2_(^^_2fi+,r)2'' J' 

(123) 

r2{x,,X2) = e-'^^^^^ ( ffc,f„, cos(Aix) - , KE,-2U)\.\/2-rW2\ 

^ ' ' 27r V 4:Al - {El - 2Q + iVy J' 

(124) 

and 

rt{xi,X2) 

~ 2n V +^/c,Ve 4A2-(i?i-21] + 2r)2'' 

(125) 

where t^^, are the two- mode single photon transmission amplitudes, and f^^, fp^ the 
two-mode single photon reflection amplitudes [Eq. fl52|) ]. Note the locations of a; and Xc in 
rt(a;i,X2) compared with ^2(2:1, 2:2) and r2{xi,X2)- 

^2(3^1, X2), T2(xi,X2), and rt(xi,a;2) represent two-photon wavefunctions in parts of the 
out-state, in which either both photons are transmitted or reflected, or one photon is trans- 
mitted while the other reflected. Experimentally, at least in principle, the magnitude of these 
wavefunctions can be measured in the setup shown in Fig. [9l where a beam splitter with a 
single-photon counter on each of the output arm, is placed at the entrance and the exit of 
the one-dimensional waveguide. In the forward (backward) direction, these photon counters 
are labeled Di, D2 (-D3, -D4), and are placed at a distance Xi, X2 from the beam splitter, 
respectively. The experiments can be carried out by injecting a weak classical beam such 
that the average numbers photons per pulse is far smaller than 2, and such that the pulse 
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repetition rate is much smaller than the inverse of the spontaneous emission lifetime. It can 
also be carried out with two-photon sources. \t2{xi,X2)\'^ corresponds to those events where 
both Di and D2 click simultaneously. The dependency on xi, X2 can be measured by varying 
the distance of the photo-detectors from the beam splitters, since \t2{xi,X2)\'^ depends only 
upon Xi — X2. Similar coincidence detection can be used in the backward direction to detect 
|r2(xi,X2)p. For |rt(xi, X2)P, one could measure the coincidence rate for Di in the forward 
direction and D4 in the backward directions. Alternatively, one could employ the Hanbury 
Brown and Twiss arrangement wherein the two photo-detectors on each side are kept at the 
same distance from the beam splitter. In this setup one measures the delay time r, which 
is proportional to Xi — X2, between two consecutive clicks on the two detectors. 

In Fig. dni we plot \t2{xi,X2)\'^, ^2(3^1, 3^2)^5 and |rt(xi,a;2)P for various total energy 
detuning 6E = E ~ 2Q, and energy difference A = {k — p)/2. Before going into details, 
we mention some general properties of \t2{xi,X2)\'^, \r2{xi,X2)\'^, and \rt{xi,X2)\'^, from the 
analytic expressions [Eqs. ( 1123p . ( ]124p . ( 1125^ ]. First of all, all \t2{xi,X2)\'^, \r2{xi,X2)\'^, and 
\rt{xi,X2)\'^ are even functions of E — 2Q and of A, thus it suffices to investigate only, say, 
the range where E — 2Q < 0, and A < 0. Also, when X2 = xi {i.e., x = 0), r2(xi,X2) is 
always zero for all E and A, i.e., the two photons are always anti-hunching in the backward 
direction. Finally, when Ai = 0, we always have \t2{ Xi, X2 — X\ )\ = V2/271, re gardless of 
the photon pair energy, E. 

We now discuss the effects of varying both SE and A. When the two incident photons 
are degenerate and on resonance with the atom, i.e., 6E = A = 0, the out-wavefucntions 
are 

t2(xi,X2) = |^e+^^^- (-e-^N/^), 
r2(a;i,a;2) = |5e-2^- (l _ e'^N/^) , 

rt{xi,X2) = ^e*2n(xi-x2)/2g-r|xi+x2|/2 ^ (_2)^ (126) 

27C 

as plotted in Fig. [TOlfa). |t2(a^i, 2:2) P decays exponentially as |a;| = \xi — X2\ becomes large, 
and thus the two transmitted photons are in a bound state. Moreover, when |x| is small, 
|t2p oc 1 — r|x| shows a cusp at x = 0, while |r2p oc does not. This should manifest in 
the measurement of the g^'^^r) function in each case. 

When the photon-pair energy is kept on resonance with the quantum impurity {5E = 0) 
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while the energy difference between the two photons, |A|, is gradually increased from zero 
to r/2, as shown in Fig. [10] (a) - (d), the peak at a; = in \t2{xi,X2)\^ reduces from its 
maximum to zero. The transmitted photons thus change from bunching to anti-bunching. 
Hence the quantum impurity can induce either an effective repulsion or attraction between 
two photons. |r2(xi,X2)P is always zero when previously mentioned. Both 

1^2(3^1, 3^2) P and |r2(xi, X2)\^ are even functions of Xi—X2- On the other hand, |rt(xi, X2)P can 
show asymmetry as a function of X1+X2 when A 7^ 0. A symmetric peak at X1+X2 = occurs 
when SE = A = (Fig. [TOlfa)). and becomes asymmetric when |A| increases (Fig. [TOT b) 
- (d)). When A 7^ 0, the maximum of |rt(xi,X2)P always occurs at xi + X2 < 0, which 
indicates the reflected photon leaves the impurity earlier than the transmitted photon. In 
addition, at 6E = 0, A 7^ 0, all the two-photon out- wavef unctions show oscillations for large 
Xi — X2 or xi + X2- On the other hand, when 6E 7^ 0, but A = 0, the oscillations at large 
Xi — X2 or xi + X2 disappear, as shown in Fig. [10] (e) and (f). 

The anti-bunching in r2{xi,X2) at x = for all 6E and A in fact has similar physical 
origin as the anti-bunching experimentally observed in resonance fluorescence from a sin- 
gle trapped ion^i. Since r2(xi,X2) arises entirely from the emission of the atom with not 
contribution from the incident light, r2{x = 0) = simply indicates that two photons can 
not be simultaneously emitted by a single atom. As a further validation of this argument, 
as well as a somewhat indirect experimental support of our theory, we note that our cal- 
culated \r2{xi,X2)\'^, as shown in Fig. [TO](a). (e), and (f), in fact agrees excellently, after 
normalization, with the experimentally measured g^'^^r) for a single trapped ion subject to 
a weak beant^. On the other hand, the predictions here for t2{xi, X2) and rt{xi, X2) involves 
interference between the incident and emitted photons and therefore represent new physical 
effects. 

The momentum distributions in each case can also be computed directly. In the forward 
direction, the momentum distribution is (again, please refer to Appendix [E] for details) 

RR{Sk2,P2\^\Skr,pi)RR = hjpi [5{ki - k2)5{pi -P2) + 5{ki -p2)5{pi - k2)] + ^B6{Ei - E2), 

(127) 

and the momentum distribution in the backward direction: 

LL{Sk2,p2\S\Ski,p^)RR = Tk^fp^ [5{ki + k2)5{pi +P2) + 5{ki +P2)S{pi + ^2)] + ^B5{Ei - E2). 

(128) 
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Define 

\k^,p^^)RL = j rfxidx2^e^'=^^^+^^'^^^4(xi)c[(a;2)|0,-), (129) 
the momentum distribution in the RL subspace is 

HL(fcf ,P2 |S|5'fei,pi)ijR = 4irpi5(A;2-^i)5(p2+Pi)+rfcitpi5(A;2-pi)5(p2 + /^i) + ^55(-£^i--£^2), 

(130) 

where 

16ir2 E1-2VL + iV 

~ Tx [4A2 - {El -2n + iVf] [4A2 - {E^ - 2^] + iVf] ' ^ 

In each momentum distribution of Eq. fll27p . (11281) . and fll30p . the delta function terms 
correspond to the uncorrelated part of the two-photon transport. The \B5{Ei — E2) term, 
however, is the signature of the strong correlation between the two photons and represents 
the background fluorescence. Specifically, B is the momentum distribution of the two pho- 
tons scattered out of the original values ki and pi. This term originates from the ee subspace, 
and gives the same contributions in the RR, LL, and RL subspaces. Fig. [11] plots normal- 
ized \B{E, Ai, A2)P as a function of Ai and A2 for various photon-pair energy E. Since the 
locations of the poles in B are at ki 2 = Pi,2 = ^ — which correspond approximately to 
either one of the photons having an energy at Q, one can picture the background fluorescence 
as one photon inelastically scattering off a composite transient object formed by the atom 
absorbing the other photon. 

The out-state wavef unctions, t2{xi,X2), r2(xi,X2), and rt{xi,X2) [Eqs. (11231) . (11241) . 
and (11251) ]. together with the corresponding momentum distributions [Eqs. (I127p . (11281) . 
(11301) ] provide a complete full quantum-mechanical description for the two-photon in-state 
\Ski,pi)RFt. scattering off a two-level system. In a classic paper, B. R. MoUow investigated 
the power spectrum of light scattered by two-level systems in a three-dimensional system, 
using a semiclassical treatment, wherein the two- level atom is driven near resonance by a 
monochromatic classical electric fielcP^. We note that in MoUow's paper, the power spec- 
trum of the scattered field, in the limit of very low incident field intensity, has exactly the 
same lineshape as the momentum distribution LL{Sk2,p2\^\Ski,pi) rb. [Eq- (]130p ] in the present 
work^''. In particular, the inelastic part of the power spectrum in Mollow's paper corresponds 
directly to the background fluorescence, In his case, however, the strength of inelastic 
scattering vanishes in the weak field limit, while in our case, strong inelastic scattering occurs 
even with only two incident photons. Therefore, the strong interference in one-dimension 
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greatly enhances the inelastic components. Also, the full quantum-mechanical treatment 
gives the correct g^'^\r) correlation function, and points out the connection between g^'^^r) 
correlation function and the background fluorescence, which could not be obtained in the 
semiclassical treatment. 

VIII. THREE PHOTON CASE 

The above procedures can be generalized to multi-photon case. For example, when there 
are three photons and one two-level system in the one-dimensional waveguide, the self- 
consistency condition becomes 



which is the generalization of Eq. (1671) . Here e{xi,X2) is the two-photon amplitude when 
one photon is absorbed and the two-level system is in the excited state. By equating the 
components, Eq. f ll32p gives a set of six self-consistent equations. These equations are exactly 
the Yang-Baxter equationsi^i^, and are connected with the integrability of the Hamiltonian. 
The details will be presented elsewhere^. 

IX. CONCLUSION AND OUTLOOK 

In this paper, we present an exact and complete solution of the transport properties 
of two-photons interacting with a single two-level system, when the photons are confined 
to a one-dimensional waveguide. Because the two-level system, at a given time, can only 
absorb one photon, the solution exhibits rich features, including, for example, the effects of 
background fluorescence and two-photon bound states. These results could be of relevance 
for many on-going quantum optics experiments. 

Also, from a formalism point of view, here we outline a general approach, based upon 
the Bethe ansatz, to solve for the transport properties of multi-particle states in a class of 
quantum-impurity problem in one- dimension. In particular, we introduce a rigorous program 
to extract the "in-" and "out-" states from the eigenstates of the interacting Hamiltonian, 
as well as a systematic approach to construct the complete scattering matrix of the system 
based upon these in- and out-states. This approach should be of general importance for a 
wide range of theoretical problems both in quantum optics and in condensed matter physics. 




(132) 
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A key observation from our solution is that the in- and out-states, as obtained from the 
standard Bethe ansatz solution, is in fact not complete, at least for the photon-Hamitonian 
(as well as the Anderson Hamiltonian in the infinite-f/ limit). The completeness of the Bethe 
ansatz solution for the interacting Hamiltonian was a subject of debate^^i^i^ in the first few 
years since the publication of the pioneering papers by Wiegmann et alM^^^^ and by Andrei 
et alA Here we note, that the proof of completeness by Schula^^, as cited by a comprehensive 
review article in this area^^, in fact only proves that completeness of the Sk,p{xi,X2) and 
Ak^p{xi,X2) states. Since the completeness of in- and out-states is particularly important 
when constructing the full scattering matrix, one needs to carefully re-examine the recent 
works of applying Bethe ansatz to the interacting resonance level system for open systems^^, 
where the crucial property of completeness of the solution is not explicitly checked. 
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APPENDIX A: COMPUTING THE IN-STATE \i) AND THE OUT-STATE \fi) 
FROM THE INTERACTING EIGENSTATE 

In this appendix, we detail the derivations of the in-state and out-state from the eigen- 
states of the interacting Hamiltonian for both the one-photon and the two-photon cases, 
which are mentioned previously in Sec. |Vl and Sec. IVI Bl respectively. Since the discussions 
are in "e" and "ee" subspaces, in this appendix, we suppress the label "e" and "ee" when 
there is no confusion. 

1. One-Photon Case 

Here, we seek to prove the forms of one-photon in-state and out-state, Eq. fH2]) . starting 
from the eigenstate \k~^)e in Eq. fHOj) . To do so, we first note that the real-space repre- 
sentation of the advanced Green's function for = J dxcKx) {—i-§^ Ce(x) and E = k, 
is 

{x, -\Gi\x', -) = (x, -I , . \x', -) = e{x' - x){+i)e^'^^--'\ (Al) 
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where 

|x,-)^ct(a;)|0,-) (A2) 

Let the out-state 1/^) be 

|/fc)= [ rfx0/(x)ct(x)|0,-), (A3) 



where (pf{x) = {x,—\fk) is the one-photon wavefunction. Projecting the Lippmann- 
Schwinger equation, Eq. (jlj) to {x,—\, we have 

= - {x, -\G^H;n,\k+). (A4) 

Inserting complete sets before and after -ffint, since Gq does not excite the impurity, and 
ifint connects {x', — | to the state |0, -|-), and does not vanish only when x' = 0, we have 



<Pf{x) = <P{x) - j dx' {x, -\G^\x', -){x', -|ifint|0, +)(+, 0|A;+) 
= (P{x)-{x,-\Gt\0,-)Vek 
= 0(x) - (^(-x)(+i)e'^"V) 



1 V 



e{-x)^= + e{x)tk^= - o{-x 



{e{-x) + ^(a;)) 4-= = 4-=. (A5) 
v27r v27r 



Thus the out-state one-photon wavefunction (pj^x) = tk^-^ = tk{x\k) for all x. In the above 



/27r 

derivations, we have used 



(x', -|i/int|0, +) = (x', -I j dxV5{x) {c\x)a + c{x)a+) |0, +) 

dx1/5(a;)(0, -|c(a;')c"^(a:)(T|0, +) 
dxV5{x)5{x' — x) 
V5{x'), (A6) 



and 



(0, = (0, +\\^J dx (f){x)c\x) + eka+j |0, -) 

= Ck. (A7) 
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Similarly, by using Eq. ([3]) and the retarded Green's function 

(x, -\G^\x\ -) = (x, -I , -) = 9ix - x'){-i)^'^^--'\ (A8) 

the in-state can be shown to be 

\k) = j dx(pi{x)c^{x)\^,-), (A9) 
with = {x\k) = for all x. 



2. Two-Photon Case 



The one-mode two-photon eigenstate l^"*") for the Hmiltonian He is computed in Sec. IVIBI 
and IVI D\ and has the following form: 

= j dxidx2g{xi,X2)^c\xi)c\x2)\%,-) + j dxe{x)c\x)\%,+), (AlO) 

with g{xi, X2) and e(x) being from either the extended Wiegmann- Andrei state in Sec. IVI Bt 
as summarized in Fig. ([7]), or the bound state in Sec. IVI Dl as summarized Fig. ([8]). The 
aim here is to prove the forms of the in- and out-state, i.e., Eq. flHHl) from Eq. flAlOp . 
Let the out-state be |/) with the following form: 

I/) = j dxidx2gf{xi,X2)^c\xi)c\x2)\%,-) (All) 

It is easy to see that 

where ^ labels the atomic state, is the identity operator in the two-photon subspace. 

Projecting (xi,X2, — | from the left to Eq. (jlj), and inserting the identity operator in the 
form of Eq. flA12p between G\ and ifjnt, we have 

gf{xi,x2) = g{xi,x2) 

— J dx'idx^dx" {xi,X2, —\G\\x'i, x'2, —){x'i,x'2, —\H\atW , +){x" , (A13) 

where we have used the fact that G\ connects ( — | only to |— ). Only the matrix element 
{x'i,x'2, — \Hi^t\x" , +) appears because of the form of i^mt- 
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We compute each matrix element in the integral. First we have 
(a;",+|i+)= / dxe{x){^,+\c{x")c\x)\^,+) 



e{x"), (AM) 



and 



~(0, -\c{x'^)c{x\) j dxV5{x) {c{x)a+ + c\x)a-) c\x")\^, +) 
V 

v_ 

-—j= {5{x\)5{x'2 — x") + 5{x'i — x")5{x'2)) ■ (A15) 
v2 



dx5{x){^, -\c{x'2)c{x[)c\x)c\x")\^, -) 

dx6{x) {S{x[ — x)6{x'2 — x") + 6{x[ — x")S{x2 — x)) 



Putting back to Eq. flA13l) . we have 



dx[dx'2dx" {xi, X2, —\G^j^\x'i, X2, — ) {5{x[)5{x'2 — x") + 5{x\ — x")5{x2)) e{x") 



gf {xi,x2) = g{xi,x2) 
V 
"72, 

= g{xu X2) - ^ y dx" {{xi,X2, |0, x", -) + {x,, X2, -1^° 0, -)) e{x") 
= g{xi,X2) - V2V j rfx"(xi,X2,-|G^|0,x",-)e(x"). (A16) 
Let {xi,X2, —\G'^\x[,x'2, —) = G'^{xi,X2;x[,x'2). Since G^^ satisfies 

{Eo - - ^e) G\ = 1, (A17) 

we have 

f d d \ II 

\Eq + i— + i— iej G\{xi, X2; x[, x'2) = {xi,X2\x[, x'2). (A18) 

Inserting the identity, 

i = ^JJ dkdp\S,,,){Sk,pl (A19) 

we have 

(^Eo + i-^+i-^-ie^ G'^^{xi,X2;x[,X2) = ^ jj dkdp {xi,X2\Sk,p){Sk,p\x[,x'2), (A20) 
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which can be solved by Fourier expanding G° (a;i, X2; x'^, x'a) using Sk,p{xi,X2): 



G\{xi,X2]x\,x'2) = {xi, X2, -\G^a\x[, X2, -) 

= jj {xi,X2\Sk,p){Sk,p\GWx\,x'2) dkdp 

= jj {xi,X2\Sk,p)Gk,p{x[,X2) dkdp. 



(A21) 



We have suppressed the idle atomic degree of freedom "— ". Inserting Eq. (]A2ip into 
Eq. flA20D . we then have 



{Eq -{k + p) - ie) Gk,p{x[, x'^ = -{Sk,p\x[,x'2) 



Gk,p{xii X2) — 



2Eo-{k + p)-te 



. {Sk,p\Xi, X2) 



(A22) 
(A23) 



■G"a{xi,X2;x[,x'2) = jj _ i^j^^^p^ _ .^ {SkJx[,X2){xi,X2\Sk,p)dkdp. (A24) 



Therefore, 
G'^(xi,X2; 0,x") 



2Eo- {k + p)-ie 
1 1 



{Sk,p\0, x") {xi,X2\Sk,p)dkdp 
2 



dkdp 



2Eo-(k+p)-ie \2'kV2 



-tpx" _j_ ^—ikx" I ^^ikxi+ipx2 _j_ ^ipxi+ikx2^ 



dkdp - 



1 



,— e-*^^"/2cosAa;V^^=cosAx 
2Eo-{k + p)-ie \ 27r 



(27r)2 
1 



dA cos Ax" cos Ax 



dE- 

^ E - Eo + ie 



JE{xc~x"/2) 



[vr {6{x + x") + 6{x - x"))] e{x"/2 - x,){-2m)e'''^^"" 



= -e{x"/2 - Xe)e*^o("=-""/') [5{x + x") + 5{x - x")] , 
where again x = Xi — X2, Xc = (xi + X2)/2. In this calculation, we have used 



dA cos Ax" cos Ax 



(A25) 



2 / dA cos Ax" cos Ax 
Jo 

POO 

/ dA (cos A(x + x") + cos A(x — x")) 
Jo 

vr (5(x + x") + 6{x — x")) . 



(A26) 
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Putting back to Eq. (lA16p . we have 

gf{xi, X2) = gixi, X2) - [e{-x/2 - Xe)e*^o("^+"/')e(-x) + 9{x/2 - x,)e'^°("^-"/')e(x)] 

V2 

= g{xi, X2) - -^V [^(-Xi)e*^o"ie(x2 - Xi) + 0(-X2)e'^°"2e(xi - X2)] . (A27) 
v2 



Using Eq. (lA27l) . one can compute the out-state directly. Here we perform the exphcit 



check for the case when is the extended Wiegmann-Andrei state in Sec. IVI CI The 
photon wavefunction g{xi,X2), and e(x), are shown in Fig. [71 The calculations are done 
separately in the four quadrants of the X1-X2 plane. 

1. In xi > and 0:2 > region, due to the step functions, 9{—Xi) and 6'(— X2) in Eq. (]A27p . 

gf{xi,X2) = g{xi,X2) = t^tpiB^e'^^^+'P''^ + A3e^?'^'i+^^^2)_ 

2. In xi < and X2 > region, since X2 — xi > 0, we have 
gfixi, X2) = fi'(xi, X2) - -j=Ve''^''''^e{x2 - xi) 

-[tpBse +hAse ) lie \k - n + iV p - n + iV /2 ) 

(i l^^pikxi+ipx2 +, A^pipxi+ikx2\ 

= t,B3.— (1 - + *.^3e"— (1 - ^_;^,r/2 ) 

= tktp (Ege^'^^^+^P^^ + A3e^P^i+*'=^2) , (A28) 
where we have used e'-^^^^ie^P^^^-xi) ^ ^i(k+p)xr^ipix2-x^) ^ ^ikxi+ipx2^ 

^iEqxi ^ik(x2—xi) _ ^i{k+p)xi^ik{x2—xi) _ ^ipxi+ikx2 
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3. In xi < X2 < region, we have 



= g{xi,X2) - -^V [e'^°"^e(x2 - xi) + e'^°"2e(xi - X2)] 
V2 



ikx\+ipx2 



ipxi+ikx2 



) 



-iT 



if RoP*P(^2-3:^i) f A „ik{x2-xi) 
^lEoxi I i^P^Sti ^ l^k^Sd \ _^ ^iEoX2 



k-n + iT/2 p-n + iT/2 

j^^^ikxi+ipx2 _j_ j^^^ipxi+ikx2^ 

t ^gg'^^^l+'P^z f^A^^ipxi+ikx2 J^^^ikxi+ipx2) j^^^ipxi+ikx2) 

k-n + iT/2 ^ p~n + iT/2 ^ p-n + iT/2 ^ k-n + iT/2 



j^^^ik{xi-X2) A^Qip{xi-X2) 



p-n + iT/2 k-Vt + iT/2) _ 



_ j^^^ikxi+ipx2 I ]^ 



iTU 



iT 



k-n + iT/2 p~n + iT/2 

iTtk iT 
p-n + iT/2 ' k-n + iT/2 



(A29) 



The xi > X2 region is obtained by gf{x2,xi) = +gf {xi,X2). 

Thus, we exphcitly demonstrate that for the interacting eigenstate of the extended 
Wiegmann-Andrei form in Sec. IVI CI the out-state photon wavefunction gf{xi,X2) is 
tfctp(i?3e*'^^^"'"*^^^ + ^3e«p^i+«'=^2')^ ^]-^g entire X1-X2 plane. This is consistent with the usual 
"read-off" of the out-state by taking the xi,X2 > region of the interacting eigenstate and 
extend to the entire X1-X2 plane. 

The in-state can be computed in exactly the same fashion by starting from Eq. ([3]) and 
shown to be gi{xi,X2) = B^e^^^^^'^^^'^ + j^^^ipxi+ikx2 ^ entire X1-X2 plane. Again, it 

is consistent with the usual read-off of the in-state by taking the xi,X2 < region of the 
interacting eigenstate and extend to the entire X1-X2 plane. Similar computations have been 
done for the two-photon bound state \Be), and the same conclusion has been reached. 



APPENDIX B: OVERLAPS OF VARIOUS STATES 



In this appendix, we summarize the properties of the two complete sets {\Sk,p)ee : k < p} 
and {\Ak,p)ee : k < p} defined in Sec. |VT1 These properties are used to normalized the 
Wiegmann-Andrei state in Sec. IVI CI as well as in the completeness check in Appendix O 
In this section, we suppress the ee label since there is no confusion. 
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We first mention the following identities: 

poo 

I COS kx dk — 7:5 (x) 
Jo 

roo y 

\ sin kx dk — V— 

Jo X 

where V denotes the Cauchy principal value. Recall that 



(Bla) 
(Bib) 



{xi,X2\Sk,p) 



1 1 

1 1 
2^7^ 



a/2 



v/2i 



{xi,X2\Ak,p) = —^sgn{x) (e^'^^ig^f^^ _ gifcx2gwi^ _ ^ sgn(a;) e^^"^ sin (Ax) 



27r 



where A = {k-p)/2^k- E/2. 

The overlap between \Ski,pi) and |Afc2,P2) is 

2 



—J V2(V2i) y dxce^^^-^'^""^ J (ixsgn(x) cos Aixsin A2X 

/•oo 

X (27r5(£^2 - El)) X 2 / dxcos Aixsin A2a; 



27r2 



= - X 6{E, 



TT 



/■c 

- ^1) X 2 / 



da; 



^ sin[(A2 + Ai)a;] + ^ sin[(A2 - Ai)a;] 



-5{E2 - El) X 



TT 

i 

TT 



A2 + A1 A2-A1 
1 



= -5(£;2 - £;i) X (2A2) ^■^2_^2- 



The overlap between and \Sk2,p2) is 



2 POO roo 

= (— )^(\/2)M dxee(^2-Ei)a;c / dxcosAixcosA2X 
27r i-00 i-00 



1 /"^ 
— X {27r5{E2 - El)) x2 / 

io 



27r2 
1 



^ cos[(Ai - A2)a;] + ^ cos[(Ai + Aa)^] 



= -5(^2 - ^1) X TT (5(Ai - A2) + S{Ai + A2)) 
TT 

= 5{E2 - El) [5{Ai - A2) + 5{Ai + A2)] 
= 5{ki - k2)6{pi -P2) + 5{ki -P2)5{k2 - Pi) 
— direct term + exchange term. 



(B2) 



(B3) 



(B4) 
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The overlap between lAfc^^pj) and |^fc2,P2) is 

(^fcl,pl|^fc2,P2) 

2^ noo POO 

= {—f{-V2i){y/2i) I dxce^^^-^'^""^ cixsin Aixsin 

■^"^ J -oo J -oo 



-oo 

r>00 



^ X {27r6{E2 - El)) X 2 dx 
-5{E2 -Ei)x7r (5(Ai - A2) - (5(Ai + A2)) 



1 1 

- cos[(Ai - A2)x] - - cos[(Ai + A2)x] 



27r2 
1 

= — ( 

TC 

= 6{E2 - El) [(5(Ai - A2) - 5(Ai + A2)] 

= 6{h - k2)6{pi - P2) - 6{h - P2)S{k2 - pi) (B5) 
= direct term — exchange term. 

Various calculations in this paper involve evaluation of overlap with the state [Eq. fl^Bl) ] 

\W,,,) = ^^^^ ^ (2A|5fe,,) + tT\A,,,)) . (B6) 

For example, 



(W^fc,p|5'fci,pi) - -j==== {2A{Sk,p\Skj^^pi) - ir{Ak^p\Sk^^pi)) 



(^6{A - Ai) + 5(A + AO - Ip^^l^^ S{E - El), 



(B7) 



2A 

~ V4A2 + r2 

where V denotes Cauchy principal value. 

In both the completeness check, as well as in the evaluation of the S-matrix, one needs 
to calculate the product {Sk2,p2\Wk^p) {Wk,p\Sk-i^^p-^) . Using Eq. (IB7p . we have 

{Sk2,pMk,p){W,,p\S,,,p,) = j^f^^iE - Ei)5{E - i?2) X I 

[5{A - Ai) + 5{A + Ai)] [5{A - A2) + 5{A + A2)] (term 1) 
+ (-0 [5{A - Ai) + 5{A + Ai)] ^^^^ (term 2) 

+ (-0 [5{A - A2) + 5{A + A2)] V^^^ (term 3) 



^.2 



+ - ^^T^ Wl \ ■ (term 4) 



A2-A2 A2-A 



(B8) 
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We now evaluate these terms. 

term 2: including the prefactor 4^^, term 2 can be simplified as 
/ r\ 4A2 1 

(--J iAfTP^ Af^^^^ " ^'^^^^ " [^^^ " + + • ^^^^ 
term 5: including the prefactor 4^q^, term 3 can be simplified as 

r^-^— ^5(i?-i?i)5(i?-i?2)[5(A-A2) + 5(A + A2)]. (BIO) 



4A2 + r2 A^-A? 



term 4 '■ in evaluating term 4, we first note the Poincare - Bertrand formula^: 

1 1 
V- — -V- 



V^T^(v—^-V—^]+n'6{X-Y)6{X-Z), (BU) 



X-Y X -Z 

Y -Z \ X ~Y X-Z 
for three arbitrary variables, X, y, and Z. Hence 



A2-Af A2-A2 



A2 - A2 V A2 - A2 A2 - a: 



7r2 



+ ,iA MA I ['^(^ - ^1) + + ^1)] [^(^ - ^2) + 5(A + A2)] (B12) 
4 Ai A2 



The terms with 5-fucntions in Eq. (]B12p . including all prefactors in Eq. (IBSp . yield 
4A2 /r\ / 7r2 [5(A-Ai) + 5(A + Ai)][5(A- A2) + 5(A + A2): 



4A2 + r2 \tt J V4|Ai||A2 

4r2 



[6{A - Ai) + 6{A + Ai)] [(5(A - A2) + 6{A + A2)] , (B13) 



4A2 + r2 

which can be combined together with term 1 in Eq. fIBSP to yield 



5{E - E^)5{E - E2) X {[5(A - Ai) + 5(A + Ai)] [(5(A - A2) + 5(A + A2)]} . (B14) 
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Therefore, the end result is 



{Sk2,P2 1 ^k,p) (M4,p| 5'fci,pi ) 
= 5{E - Ei)6{E - E2) X 

[S{A - Ai) + 5{A + Ai)] [S{A - A2) + S{A + A2)] 

+ (-D4ft^A?^t^(^-^^)^^(^ + ^^)] 
+ B) 4^ ^ Af^ - + + 

/r\^ 4A2 1/1 1 



APPENDIX C: COMPLETENESS CHECK 

In this appendix, we carry out the exphcit check of the completeness of the eigenstates 
{|Wfc,p), \Be)} in Sec. IVI CI Again, since the discussions below are in the ee subspace, we 
omit the subscript when there is no confusion. 

As noted in Sec. IVI C| to check whether {|W^fc,p) : k < p} is complete, one could start with 
an arbitrary state, for example, l^'^^^p^), project out all |Vrfe,p) components and calculate 

li^fci.pi) = \Ski,pi) —''^^{Wk,p\Ski,p-^)\Wk,p). (CI) 

k<p 

If the set {iVTfc^p) : k < p} were complete, such a computation should yield l^^j^p-^) = for 
arbitrary ki and pi. To calculate |5fci,pi), we first project |5fci,pi) to {Sk2,p2\- 

{^k2,P2\^ki,pi) = {Sk2,p2\^ki,pi) ~''^^{^k,p\Ski,pi){Sk2,p2\'^k,p)- (C2) 

k<p 

The first term in the right hand side is [from Eq. (lB4p ] 

{Sk2,p2\Ski,p,) = 5{ki - k2)5{pi -P2) + 5{ki -^2)^(^1 - k2), (C3) 

while in the second term, the restriction k < p can be dropped, using the symmetry property 
of \Wk,p), Eq. iib!): 

{Sk2,P2\^k,p){Wk,p\Ski,pi) = -''^{Sk2,p2\Wk,p){Wk,p\SkT^^p^). (C4) 
k,p,k<p k,p 
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Using Eq. (IBISP for {Sk2,p2\Wk,p) (Wk^plSkj^^p-^^) , the second term becomes 

9 y^A^k2,P2\Wk,p){Wk,p\Sk^^p^) 



k,p 



dE 



dA6{E - Ei)6{E - E2) x 



[6{A - Ai) + 6{A + Ai)] [6{A - A2) + 6{A + A2)] 



+ 



r 



AAj 



V- 



1 



TT y 4A2 + r2 A2 - A 



[5(A- Ai)+5(A + Ai)] 



^ r 



4A2 



1 



TT ; 4A2 + r2 A2 - A? 



r 



4A2 



+ 



The first term in Eq. fICSp yields 



[5(A- A2) +5(A + A2)] 
1 /„ 1 „ 1 



vr y 4A2 + r2 Al- Al\ A2 - A? A2 - A 



V- 



-V- 



(C5) 



- E2) [5(Ai - A2) + 5(Ai + A2)] 

- {Sk2,p2\Skx,p^ ■ 



(C6) 



The integrations of the second and third terms are straightforward, and the sum of both 
terms give 



4r- 



1 



;5{E^-E2). 



TT 4A2 + r24A2 + r2"^^' ^'^^^ 

The last term can be calculated using a contour integral. The only non-vanishing contri- 
bution comes from the pole at A = +ir/2, when the integration contour is chosen to be 
completed in the upper half plane. Hence the integration yields 

2 / -n2 



1 V 
-2m - 

2 y-n 

4r3 



V- 



1 



1 



TT 4A2 + r2 4A2 + r2 

The final result therefore is 



4«r ; A2 - A2 V(^r/2)2 - ai (ir/2)2 - a2 



5(^1 - E2 



(C8) 



{Sk2,P2\^k,p)(Wk,p\Ski,p^) 

k,p,k<p 

5{ki - k2)S{pi -P2) + S{ki -p2)S{k2 -pi) 

Q-p3\ 1 1 

+ I 1 ... . ... . ^J iE, - E2). 



TT y 4A2 + r2 4A2 + r2 



(C9) 
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Thus 

^— ) a. ^T-J^n. - ^2). (CIO) 

Since (5'fc2,p2l^fei,pi) 7^ 0, this directly proves that the set : /c < p} is incomplete. 

A very important observation regarding Eq. (IClOp is that, independent of the choice of 
fci, pi, the resulting state {|5fci,pi)} calculated in Eq. (IClOP is always proportional to the 
same state 

2 /'+00 r-\-oo 

l^kupi) = 9 / ^-^2 / '^A2(5'fe2,p2|(5fci,pi)|S'fc2,p2) 

/H-oo /»+oo -j^ 

= (Cll) 

Therefore, the set {|5fci,pi)} in fact forms a one-dimensional Hilbert space. Thus, only one 
extra state \Be) is needed in order to span the two-photon Hilbert space. Since 

= - E'), (C12) 

this extra state, when normalized, is 



This concludes the proof that {|W^fc,p), \Be)} forms a complete basis of the two-photon 
Hilbert space. 

To see the physical meaning of \Be), we rewrite Eq. (]C13|) as 



\Be) = j dxcdxBE{xuX2)^c\xi)c\x^)\(D,~), (C14) 

with 

/r 

BEix^,X2) = Be{x,,x) = e^^-'^4=e-^l^l. (C15) 

v47r 

In the above derivation, we have used 

Thus, the state \Be) in fact defines a two-photon bound state. 
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APPENDIX D: DERIVATIONS OF S-MATRIX 



In this appendix, we provide the detailed calculations of Eq. flllSp . the matrix element 
of the S-matrix in the ee subspace, See- Again, since the discussions below are in the ee 
subspace, we omit the subscript when there is no confusion. 

The S-matrix in the ee subspace. See is defined in Eq. (11171) : 

See = J2^ktp\W,,p){W,,p\ +Y,tE\BE){BE\. (Dl) 

k<p E 

Our goal is to compute the matrix element, Eq. (11181) . 

{Sk2,P2\^ee\Ski,p^) = {Sk2,P2\ ^^^iktplWk^p) {Wk,p\^ \Ski,pi) + {Sk2,p2\ ^^^iE\BE){BE\^ \Ski,pi)- 

(D2) 

We will compute the two terms on the right hand side separately in the following. 
Dl. First term of Eq. (|D2]) 

Using Eq. flBlSp for (5'fc2,p2 l^fc,p) (^fc.pl'S'fci.pi), the first term in Eq. flD2p becomes 

tktp{Sk2,p2\Wk,p) {Wk^p\Ski,pi) 



k<p 

1 

2 



9 tk'tp{Sk2,p2\'^k,p) (W4,p|5'a:i,pi) 



k,p 



oo r+oo 

dE / dA 5{E - Ei)6{E - E2) x tktp 

—00 J ~oo 



\S{^ - Ai) + ^(A + Ai)] |<5(A - A2) + <5(A + Aj)] 

DiA!TF^AFA!l*<^-^'' + *<^ + ^')l 
+ I^)'4A^^aFa|(^a^-^aJ^I^- 
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Note that 



tktp — tA+E/2t-A+E/2 

'A + {E- 2VL)/2 - iT/2\ f-A + {E- 2VL)/2 - iV/2 



A + {E- 2n) /2 + iT/2j \-A + {E - 2fi) /2 + iV /2 
A + {E- 2VL)/2 - iT/2\ [A-{E- 2VL)/2 + iT/2 



A + {E- 2VL)/2 + iV/2 J \A - {E - 2VL) /2 - iV/2 
the first term in Eq. fIDSP can be evaluated as 



(D4) 



tk.t^MEi - E2) [5(Ai - A2) + 5(Ai - A2)] 
= tkjpj{ki - k2)5ipi -P2) +tk^tp^5{ki -P2)5{pi - k2). (D5) 

The second and the third term in Eq. (lD3p can be combined to give 
r / 4A^ 1 4A^ 1 \ 

In evaluating the fourth term in Eq. (]D3p . we use a contour integral, by completing the 
contour in the upper half plane. The only poles that give non-zero contributions are located 
at A = - 2^ + zr)/2, and A = zr/2. The result is 

Ei-2VL~ 2iT 1 1 



5{Ei -E2) y E1-2Q + 2tr 4AI + r2 4A2 + r2 
ler^ Ei-2Q + iT 1 1 

~ Ei-2n + 2iT 4Af - {Ei - 2^] + iV)^ AAj - (^1 - 2fi + iVy 

Summing Eqs. (IDSp . (]D6p . and (]D7p together, we thus have 

{Sk2,P2\ I ^4^p|W4,p)(W^fc,p| I \Ski,pi) 
\k<p / 

= tk,tp^5{ki - k2)5{pi -P2) +tk^tp^5{ki -P2)5{pi - k2) 

r / 4A2 1 4A2 1 

+ -S{Ei - E2) tk^tp^ -no^Tl ^ + ^fci^Pi.A2 , -no^ 



(D7) 



+ - F 1 El - 2^] - 2zr i i 

^ ^ ' '^V ^ El - 2fi + 2^r 4A2 + r2 4A2 + p 

ler^ Ei-2n + iT i i 

~ Ei~2n + 2iT 4Af - (Ei -2n + iT)^ AAj - {Ei - 2^] + iVy 



(D8) 
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D2. Second term of Eq. dD2]) 



We first evaluate the overlap between \Be) and \Ski^p^): 



{BE\Ski,pi) 



/CO /'OO 
-oo J — oo 



1 1 



1 1 



_27rv/2 ^ : 



47r 27r 
^ - E2 



^27r 

4r 



iAi-r/2 -iAi-r/2 



Thus, the second term of Eq. (lD2p . the bound state contribution, is 



4r 



(^-2fi + 2^) 4AfTP 4A1TP^^^ " ^'^^^^ - 



E1-2Q- 2iT 1 



AE, - E2] 



(D9) 



(DIO) 



TT Ei-2n + 2iT 4Af + 12 4Ai + r2 
Summing Eqs. fIDSp and fIDlOp . we obtain the S-matrix {Sk2,p2\See\Ski,pi) , Eq. (11180 : 

{Sk2^P2\See\Sk^.p^) = tkJp^S{ki-k2)5{pi-p2)+tk,tp^Siki-p2)5{k2-pi)+BS{Ei-E2), (Dll) 



with 



B 



Ei-2n + iV 



TT [4Af - (^1 -2n + tvy] [4Ai - {El -2n + iv)^] 



(D12) 



APPENDIX E: DERIVATIONS OF TWO-MODE OUT-STATE 



In this appendix, we present the details of the derivations of the two-mode out-state two- 
photon wavefunciton, t2{xi,X2) [Eq. fll23p ]. r2{xi,X2) [Eq. fll24p ]. andrt(xi,X2) [Eq. fll25p ]. 
The in-state is a state of two right-going photons, \Sk-i,p^)RR. We first decompose the 
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in-state into ee, oo, and eo subspaces: 



I in) = \Skj^^p^)RR 



27rv2 v2 



+ 1 dx.dx,^^ ^^rkx.+^px. + (4(Xi)ct(x2) + 4 (xOc^X^)) |0, -) 

= l\Sk,,P.)ee + l\S,„,^)oo + 1^1 ^^^kx,+^px, ^ ^kx.Mpx,^ ^ (^^^(x^) |0, -) 

= |in)ee + |in)oo + |in)eo (El) 

Employing the decomposition relation, Eq. (!32|) . we carry out the calculations in each 
individual subspace: 

S|in)ee = See|in)ee = See-|5,„,,)ee 



ee ee^S J7;2,A2 I ^ee | 'S'fcj^^pj^ )ee 

-B2,A2<0 

- f 4i^pi|5'fci,pi)ee + ^ -B|5'£;i_A2)ee j 
\ A2<0 / 

If 1 

- y rfxi(ix2ee(a;c,a;|out)ee--^4(xi)4(x2)|0, -) 



If 1 

= -y cixi(ix2 0ee(a;i,X2)--^4(3;i)4(3;2)|0, -), (E2) 



0ee(a;i, 3:2) is the out-state wavefunction in ee subspace, Eq. fll2Up . 
S I in) 00 = Soo|in)oo = Sqo— l^fe-^ p-^)oo 



1 /" 1 

- j dXidX2oo{Xi,X2\Sk^^p^)oo-^cl{Xi)cl{x2)\%,-) 
If 1 

- J dxidx2Sk,,p,^cl{xi)cl{x2)\^,-) (E3) 
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and 



S|in)eo = Seo|in)f 



1 1 



1 



-= I dxidx2 ^ f< 

2^ J ^271^2^ 



ikxi+ipx2 ^ikx2+ipxi\ (^\{^.^\^\ 



) (4(xi)ct(x2)+ct(xi)4(x2)) |0) 



(E4) 



Using the transformation formula, Eq. (1271) . and collect terms according to the operators 
in Eq. (11221) . we obtain t2{xi,X2), T2(xi,X2), and rt{xi,X2): 

t2{Xi,X2) = ^ [(l)ee{Xi,X2) + Sk^^p^{Xi, X2) + (tfc + tp) Sk^^p^{xi, X2)] 



1 

4 

1 
4 



tkitpiSki,piiXi, X2) + BSei,A2{^cj + 'S'/fc^^pj (Xi, X2) + (tfe + tp)5'fe^^p^ (Xi, X2) 

A2<£;i/2 



,A2 (-^C) 

A2<Si/2 



A2<Bi/2 

72 /- ^ ... 



AEiXc 



— ( ^i^pi COS (Aix) - 



4A2 - (El - 21] + ir)2 



^i(£;i-2C)|3;|/2-r|x|/2 



(E5) 



(E6) 



r2(Xi,X2) = - [(/)ee(-Xi, -X2) + Sk^^p^{-Xi, -X2) - {tk + tp)Sk^^p^{-Xi, -X2)] 



1 

4 



(1 - tfcj(l - tpJS'fci,p,(-xi, -X2) + ^ 5S'ij,,A2(-a;c, -a;) 

A2<Si/2 



'"fci'^pi5'fc,,pi(-a;i, -X2) + ^ ^ SS'£i,A2(-a;c, -a:) 

A2<£i/2 



,i{Ei-2Q.)\x\/2-T\x\/2 



(E7) 
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and 



rt{xi,X2) 
1 



4^2 



(j)ee{xi, -X2) - Sk^^p^{xi, -X2) + {tp - tk) 



-J_(, 



,i(Ei-2f2)|Xc|-r|Xe| 



(E8) 



The momentum distributions can be computed directly. In the forward direction, 



RR{Sk2,p2\^\^ki,pi)RR 

{RR{Sk2,p2\^'^^) RR + RR{Sk2,P2\^ki,pi)RR + tkitp^RR{Sk2,p2\Ski,pi) Rr) 



= - [{h, + l){tp, + 1) {6{h - k2)S{p, - P2) + 6{h - p2)S{pi - k2)) + B5{E^ - E2) 
= hjp^ {S{h - k2)6{p, - P2) + 5{ki - p2)6{pi - k2)) + \b6{E^ - E2). 



(E9) 



In the backward direction: 



LL{Sk2,p2\^\Ski,pi) RR 

-- J dxidx2Sl^ p^{xi,X2)r2{xi,X2) 



= - [(1 - tk,)il - tp,) i5iki + k2)5ip, +P2) + 5iki +p2)Sipi + k2)) + B5{E^ - E2)\ 



'^k-J'pi 

{d{ki + k2)dipi + P2) + S{ki+P2)dipi + k2)) + -B5{E^ - E2), 



(ElO) 



while in the RL subspace: 



RL{k2,P2\^\Ski,pi)RR 



= J dxidx2 ^i-e^'^^^i+^f^^^^ rt{xi,X2) 



■ik,fpAk2 - ki)5{p2+pi)+fk,tp,5{k2-pi)5{p2 + /ci) + -B5{Ei - E2). (Ell) 



In the above calculations, we have adopted the following sign convention for the left-moving 
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photons: 

dx^=c{{x)\^,-), with < 0. 
v27r 



\Sk,p)LL = jj dxidx2 (e^^-i+*P-2 + e^k.,+^p.,^ -Lc[(xi)c[(x2)|0, -), with k,p<0. 

\k'',p'^)RL rfxidx2^e^'^"i+^P^24(xi)4(x2)|0, -), with A; > 0, p < 0. (E12) 
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FIG. 1: (Color online) Pictorial representation of the Hamiltonian of Eq. ([1]). The dynamics of the 
free quantum many-particles are described by Hp, which can have more than one incident particle. 
The dynamics of the impurity (the "atom") is described by Ha- The interactions between the 
quantum particles and the impurity is described by -f^int- 
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FIG. 2: The evolution of the state of the system with adiabatic switching. \i) and are free 
states governed by the free Hamiltonian, Hq. When the interaction is adiabatically switched on 
from the distant past {t — > — cxd) to its full strength i^int at t = 0, the state evolves from asymptotic 
state \J^{t — > — cxD)|i) to {i'^), which is governed by the full interacting Hamiltonian H = i^o + i^int- 
When the interaction is adiabatically switched off in the remote future {t — > +00), the state 
asymptotically approaches U°(t —>■ +oo)|/j). 
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FIG. 3: A quantum particle incident upon a delta potential barrier characterized by Vod{x). 
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FIG. 4: (Color online) Schematics of the system. A two-level system is coupled to a one-dimensional 
continuum in which the photons, shown as wiggly waves, propagate in each direction. 
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FIG. 5: The single-photon transmission and the reflection spectrum. is indicated by the black 
curve, and |ffcp is denoted by the gray curve. The full width at half maximum for |ffcp is F. 
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FIG. 6: The xi-axis, a;2-axis, and xi = X2 dissect the X1-X2 coordinate plane into six regions, labeled 
by the numbers in circle. When given g{xi,X2) in region 3 (lightly-shaded area), the boundary 
condition is imposed to obtain g{xi,X2) in other regions, as denoted by the arrows. g{xi,X2) in 
xi > X2 region (darkly-shaded area) is obtained from g{xi,X2) in X2 > xi region by the boson 
statistics. 
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rr < X > X 



B3 _k-p-ir 
A3 " k-p + iV 

FIG. 7: The wavefunction of the interacting eigenstate \i^) of He as constructed from the standard 
Bethe ansatz approach. The shaded region is obtained by symmetry, i.e., g{x2,xi) = +g{xi,X2). 
Also shown is 6(2;) for ah x. 
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e[x) 

2^/2V iEx/2-^Tx/2 2v/2F iEx/2~Tx/2 

E-2n + 2iV E-2n + 2iT 

► 

X < a: > X 

_ E-^2n- 2iT 

FIG. 8: The wavefunction of the interacting eigenstate of cosisting of two-photon bound 
state. The shaded region is obtained by symmetry, i.e., g{x2,xi) = +g{xi,X2). Also shown is e(x) 
for all X. 
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FIG. 9: (Color online). Schematic experimental setups for concurrence measurements of (a) 
\t2{xi, X2)\'^ , (b) |r"2(2;i, X2)P, and (c) |rt(xi,X2)p. Di, D2 are photo-detectors with adjustable 
positions. BS, beam splitter. The "=" symbol inside the one-dimensional waveguide denotes the 
two-level quantum impurity. 
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FIG. 10: \t2{x\,X2)^ , |''2(iCi, iC2)P, and \rt{x\^X2)\'^ for various photon-pair energy detuning 5E = 
E - 2n, and energy difference A. (a) SE = 0, A = 0, (b) SE = 0, A = -0.3T, (c) SE = 0, 
A = -0.4r, (d) SE = 0, A = -0.5r, and (e) 6E = -l.bF, A = 0, (f) SE = -2.5r, A = 0. 
X = Tx/2. 
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FIG. 11: (Color online) Background fluorescence as a function of Ai and A2 at various energy, (a) 
^ = 0. (b) ^ = 2. (c) E = A. (d) E = 6. B = {T/2)B, E = {E - 2n)/{T/2), and A = A/(r/2). 
For any given E, the in- and out-states can be completely specified by one quadrant in the A1-A2 
plane. 
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